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Abstract 



, Partition functions, also known as homomorphism functions, form a rich family 

>■• ' of graph invariants that contain combinatorial invariants such as the number of k- 

C/3 , colourings or the number of independent sets of a graph and also the partition functions 

of certain "spin glass" models of statistical physics such as the Ising model. 

Building on earlier work by Dyer and Greenhill 9J and Bulatov and Grohe [6], 
CO ' we completely classify the computational complexity of partition functions. Our main 

^ , result is a dichotomy theorem stating that every partition function is either computable 

in polynomial time or #P-complete. Partition functions are described by symmetric 
Qs , matrices with real entries, and we prove that it is decidable in polynomial time in terms 

of the matrix whether a given partition function is in polynomial time or #P-complete. 
While in general it is very complicated to give an explicit algebraic or combinatorial 
f~^ ' description of the tractable cases, for partition functions described by a Hadamard 

00 , matrices — these turn out to be central in our proofs — we obtain a simple algebraic 

^D ■ tractability criterion, which says that the tractable cases are those "representable" by 

a quadratic polynomial over the field F2. 

X 



* Partly funded by the EPSRC grant "The complexity of counting in constraint satisfaction problems". 

and by the Deutsche Forschungsgemeinschaft within the research training group 'Methods for Discrete 

Structures' (GRK 1408) 
^Department of Computer Science, University of Liverpool, Liverpool L69 3BX, UK 
"'"Institut fiir Informatik, Humboldt-Universitat zu Berlin, 10099 Berlin, Germany 
^School of Mathematical Sciences, Queen Mary, University of London, Mile End Road, London El 4NS, 

UK 
^Institut fiir Informatik, Humboldt-Universitat zu Berlin, 10099 Berlin, Germany 



1 Introduction 

We study the complexity of a family of graph invariants known as partition functions or 
homomorphism functions (see, for example, |14 1 I2 H [22]). Many natural graph invariants 
can be expressed as homomorphism functions, among them the number of /c-colourings, 
the number of independent sets, and the number of nowhere-zero A:-flows of a graph. The 
functions also appear as the partition functions of certain "spin-glass" models of statistical 
physics such as the Ising model or the g-state Potts model. 

Let A G ]^"^x™- l^e a symmetric real matrix with entries Aij. The partition function Za 
associates with every graph G = {V, E) the real number 

Za{G) = X] n Mn),i(v)- 

g:y^[m] {u,v}eE 

We refer to the row and column indices of the matrix, which are elements of [m] := 
{1, . . . ,m}, as spins. We use the term configuration to refer to a mapping S, '■ V ^ [m] 
assigning a spin to each vertex of the graph. To avoid difficulties with models of real number 
computation, throughtout this paper we restrict our attention to algebraic numbers. Let 
Ma denote the set of algebraic real numbers o 

Our main result is a dichotomy theorem stating that for every symmetric matrix A G 
j^mxm ^YiQ partition function Z^ is either computable in polynomial time or ^^P-hard. 
This extends earlier results by Dyer and Greenhill [9], who proved the dichotomy for 0-1- 
matrices, and Bulatov and Grohe 0, who proved it for nonnegative matrices. Therefore, 
in this paper we are mainly interested in matrices with negative entries. 

Examples 

In the following, let G = {V, E) be a graph with N vertices. Consider the matrices 

'0 1' 
S = ( 1 and C?, 

It is not hard to see that Zs{G) is the number of independent sets of a graph G and 
Zc:j{G) is the number of 3-colourings of G. More generally, if A is the adjacency matrix 
of a graph H then Za(G) is the number of homomorphisms from G to H. Here we allow 
H to have loops and parallel edges; the entry Aij in the adjacency matrix is the number 
of edges from vertex i to vertex j. 

Let us turn to matrices with negative entries. Consider 

H.= (\ _\). (1.1) 




'^ There is a problem with the treatment of real numbers in [B], but aU results stated in [B] are valid for 
algebraic real numbers. We use a standard representation of algebraic numbers by polynomials and 
standard Turing machines as our underlying model of computation. 



Then 2Zh2{G) + 2 is the number of induced subgraphs of G with an even number of 
edges. Hence up to a simple transformation, Zh2 counts induced subgraphs with an even 
number of edges. To see this, observe that for every configuration ^ : y ^ [2] the term 
riiu u>ge ^£.{u),^(v) is 1 if the subgraph of G induced by ^^"^(2) has an even number of edges 
and —1 otherwise. Note that H2 is the simplest nontrivial Hadamard matrix. Hadamard 
matrices will play a central role in this paper. Another simple example is the matrix 

y -1 1 

It is a nice exercise to verify that for connected G the number Zu{G) is 2 if G is Eulerian 
and otherwise. 

A less obvious example of a counting function that can be expressed in terms of a 
partition function is the number of nowhere-zero fe-flows of a graph. It can be shown that 
the number of nowhere-zero A;-flows of a graph G with N vertices is k~ ■ Zp^{G), where 
-Ffc is the k X k matrix with {k — l)s on the diagonal and —Is everywhere else. This is a 
special case of a more general connection between partition functions for matrices A with 
diagonal entries d and off diagonal entries c and certain values of the Tutte polynomial. 
This well-known connection can be derived by establishing certain contraction-deletion 
identities for the partition functions. For example, it follows from [24\ Equations (3.5.4)] 
and t23i Equation (2.26) and (2.9)] 

Complexity 

Like the complexity of graph polynomials [21 [lU [181 [20] and constraint satisfaction prob- 
lems [H [31 [H [5l [121 [151 [E]) which are both closely related to our partition functions, the 
complexity of partition functions has already received quite a bit of a attention. Dyer and 
Greenhill [9j studied the complexity of counting homomorphisms from a given graph G to 
a fixed graph H without parallel edges. (Homomorphisms from G to H are also known as 
H -colourings of G.) They proved that the problem is in polynomial time if every connected 
component of H is either a complete graph with a loop at every vertex or a complete bi- 
partite graph, and the problem is #P-hard otherwise. Note that, in particular, this gives 
a complete classification of the complexity of computing Z^ for symmetric 0-1-matrices 
A. Bulatov and Grohe [6] extended this to symmetric nonnegative matrices. To state the 
result, it is convenient to introduce the notion of a block of a matrix A. To define the 
blocks of A, it is best to view A as the adjacency matrix of a graph with weighted edges; 
then each non-bipartite connected component of this graph corresponds to one block and 
each bipartite connected component corresponds to two blocks. A formal definition will 
be given below. Bulatov and Grohe [6] proved that computing the function Za is in poly- 
nomial time if the row rank of every block of A is 1 and ^^P-hard otherwise. The problem 
for matrices with negative entries was left open. In particular, Bulatov and Grohe asked 
for the complexity of the partition function Zh^ for the matrix H2 introduced in (jl.ip . 
Note that H2 is a matrix with one block of row rank 2. As we shall see, Zh2 is computable 
in polynomial time. Hence the complexity classification of Bulatov and Grohe does not 



extend to matrices with negative entries. Nevertheless, we obtain a dichotomy, and this is 
our main result. 

Results and outline of the proofs 

Our main theorem is the following. 

Theorem 1.1 (Dichotomy Theorem). Let A G M™^™ be a symmetric matrix. Then 
the function Za either can be computed in polynomial time or is #P-hard. 

Furthermore, there is a polynomial time algorithm that, given the matrix A, decides 
whether Za is in polynomial time or #P-hard. 

Let us call a matrix A tractable if Za can be computed in polynomial time and hard 
if computing Za is ^^P-hard. Then the Dichotomy Theorem states that every symmetric 
matrix with entries in Ma is either tractable or hard. The classification of matrices into 
tractable and hard ones can be made explicit, but is very complicated and does not give 
any real insights. Very roughly, a matrix A is tractable if each of its blocks can be written 
as a tensor product of a positive matrix of row rank 1 and a tractable Hadamard matrix. 
Unfortunately, the real classification is not that simple, but for now let us focus on tractable 
Hadamard matrices. Recall that a Hadamard matrix is a square matrix H with entries 
from { — 1, 1} such that H ■ H is a diagonal matrix. Let H ^ { — !> i|"X" '\qq a symmetric 
n X n Hadamard matrix with n = 2 . Let p : Fg ^ [n] be a bijective mapping, which we 
call an index mapping. We say that a multivariate polynomial h{Xi, . . . , X^, Yi, . . . , Y^) 
over F2 symmetrically represents H with respect to p if, for all x = (xi, . . . ,Xk),y = 
(yi, . . . , yk) G F2, it holds that 

h{xi,...,Xk,yi,...,yk) = 1 ^=^ i;rp(x),p(y) = -1. 

For example, the F2-polynomial h2{Xi,Yi) = Xi ■ Yi symmetrically represents the 
matrix H2 with respect to the index mapping p{xi) = xi + 1. The F2-polynomial 
/14 (Xi , X2 , Yi , Y2) = Xi • Y2 (B X2 • Yi symmetrically represents the matrix 



Ha 



with respect to the index mapping p{xi,X2) = 2-xi+X2 + l. The qualifier "symmetrically" 
in "symmetrically represents" indicates that the same index mapping is applied to both 
X and y. We will need to consider asymmetric representations later. Note that we can 
only represent a matrix H £ { — 1, l}"-^" by an F2-polynomial in this way if n is a power of 
2. In this case, for every index mapping p there is a unique F2-polynomial symmetrically 
representing h with respect to p. We say that H has a quadratic representation if there 
is an index mapping p and an F2-polynomial h of degree at most 2 that symmetrically 
represents H with respect to p. Our dichotomy theorem for Hadamard matrices is as 
follows. 
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Theorem 1.2 (Complexity Classification for Hadamard Matrices). A symmetric 
Hadam^ard matrix H is tractable if it has a quadratic representation and hard otherwise. 

Hence, in particular, the matrices H2 and H4^ are tractable. The tractability part of 
Theorem 11.21 is an easy consequence of the fact that counting the number of solutions of a 
quadratic equation over F2 (or any other finite field) is in polynomial time (see |131 [T9] ) . 
The difficulty in proving the hardness part is that the degree of a polynomial representing 
a Hadamard matrix is not invariant under the choice of the index mapping p. However, for 
normalised Hadamard matrices, that is, Hadamard matrices whose first row and column 
consists entirely of +ls, we can show that either they are hard or they can be written 
as an iterated tensor product of the two simple Hadamard matrices H2 and -^4. This 
gives us a canonical index mapping and hence a canonical representation by a quadratic 
F2-polynomial. Unfortunately, we could not find a direct reduction from arbitrary to 
normalised Hadamard matrices. To get a reduction, we first need to work with a gener- 
alisation of partition functions. If we view the matrix A defining a partition function as 
an edge-weighted graph, then this is the natural generalisation to graphs with edge and 
vertex weights. Let A £ M^^"* be a symmetric matrix and D G M™^"^ a diagonal matrix, 
which may be viewed as assigning the weight Z?i_j to each vertex i. We define the partition 
function Za,d by 

ZaMG) = Z] n ^€(«).€{'') ■ n ^CM.CM' 

^■.V-'[m] {u,v}&E veV 

for every graph G = {V,E). As a matter of fact, we need a further generalisation that 
takes into account that vertices of even and odd degree behave differently when it comes 
to negative edge weights. For a symmetric matrix A £ R™^™ and two diagonal matrices 
Z),OeM^^™ welet 

za,d,o{g)= yi n ^?(«).ew n ^i;{v),ii{v) ■ n '^^m.^M' 

^■.V~*[m] {u,v}eE v£V v£V 

deg(t') is even deg(i)) is odd 

for every graph G = {V,E). We call Za.d.o the parity- distinguishing partition function 
(pdpf) defined by A,D,0. We show that the problem of computing Za,d,o{G) is always 
either polynomial-time solvable or ^^P-hard, and we call a triple (^4, D, O) tractable or 
hard accordingly. Obviously, \i D = O = Im are identity matrices, then we have Za = 

Za,d = Za,d,o- 

Returning to the outline of the proof of Theorem II. 2^ we can show that, for every 
Hadamard matrix ff, either H is hard or there is a normalised Hadamard matrix H' 
and diagonal matrices D' , O' such that computing Zh is polynomial time equivalent to 
computing Zhijji^o'- Actually, it turns out that we may assume D' to be an identity 
matrix and O' to be a diagonal matrix with entries 0, 1 only. For the normalised matrix 
H' we have a canonical index mapping, and we can use this to represent the matrices 
D' and O' over F2. Then we obtain a tractability criterion that essentially says that 
(H' ,D',0') is tractable if the representation of H' is quadratic and that of O' is linear 
(remember that D' is an identity matrix, which we do not have to worry about). 



For the proof of the Dichotomy Theorem II. H we actually need an extension of Theo- 
rem [L2] that states a dichotomy for parity-distinguishing partition functions Za,d,0) where 
j4 is a "bipartisation" of a Hadamard matrix (this notion will be defined later) . The proof 
sketched above can be generalised to give this extension. Then to prove the Dichotomy 
Theorem, we first reduce the problem of computing Z^ to the problem of computing Zq 
for the connected components C of A. The next step is to eliminate duplicate rows and 
columns in the matrix, which can be done at the price of introducing vertex weights. Us- 
ing the classification theorem for nonnegative matrices and some gadgetry, from there we 
get the desired reduction to parity-distinguishing partition functions for bipartisations of 
Hadamard matrices. 

Let us finally mention that our proof shows that the Dichotomy Theorem not only holds 
for simple partition functions Z^, but also for vertex- weighted and parity-distinguishing 
partition functions. 

Preliminaries 

Let A G M™^" be an (m x n)-matrix. The entries of A are denoted by ^ij. The ith row 
of A is denoted by ^j,*, and the jth column by ^*j. By abs(A) we denote the matrix 
obtained from A by taking the absolute value of each entry in A. 

Let Im be the m x m identity matrix and let Im;A be the m x m matrix that is all zero 
except that Ijj = 1 for j S A. 

The Hadamard product C of two m x n matrices A and B, written C = A o B, is 
the m X n component-wise product in which Cij = AijBij. —A denotes the Hadamard 
product of A and the matrix in which every entry is —1. 

We write {u,v) to denote the inner product (or dot product) of two vectors in M^. 

Recall that the tensor product (or Kronecker product) of an r x s matrix B and antxu 
matrix C is an rt x su matrix B ® C. For A; G [r], i G [t], £ € [s] and j £ [u], we have 
(5(8)C)(fc_i)t_i_j^(£_i)„+j = Bk/Cij. It is sometimes useful to think of the product in terms 
of rs "blocks" or "tiles" of size t x u. 



B^C 



BiiC . . . BisC 



BrlO . . . BygO 



For index sets / C [m], J C [n], we let Ai^j be the (|/| x \J\)-suhmatrix with entries Ai^j 
for i € I, j € J. The matrix A is indecomposable if there are no index sets / C [m], J C [n] 
such that {I, J) / (0,0), (/, J) / ([m], [n]) and Aij = for all (i, j) G (([m] \ I) x j) U 
(/ X ([n] \ J)) . Note that, in particular, an indecomposable matrix has at least one nonzero 
entry. The blocks of a matrix are the maximal indecomposable submatrices. For every 
symmetric matrix A £ M"^" we can define a graph G with vertex set [n] and edge set 
{{i, j} I j4j j 7^ 0|. We call the matrix A bipartite if the graph G is bipartite. We call A 
connected if the graph G is connected. The connected components of A are the maximal 
submatrices Ac^c such that G[C], the subgraph of G induced by C C [n], is a connected 
component. If the connected component G[G] is not bipartite then Ac,c is a block of A. If 



the connected component G[C] is bipartite and contains an edge then Ac^c has the form 
( rT n ) ' '^'^6^6 i? is a block of A. Furthermore, all blocks of A arise from connected 

components in this way. 

For two Counting Problems / and g, we write f < g ii there is a polynomial time Turing 
reduction from f to g. If / < 5 and g < f holds, we write f = g. For a symmetric matrix 
A and diagonal matrices D, O of the same size, EVAL(A, D, O) (EVAL(^, D), EVAL(^)) 
denotes the problem of computing Za,d,o{G) {Za,d{G), Za{G), respectively) for an input 
graph G (which need not be a simple graph - it may have loops and/or multi-edges). 

2 Hadamard matrices 

The main focus of this section is to prove Theorem 12.21 below which is a strengthened 
version of Theorem 11.21 Suppose that H is a.n n x n Hadamard matrix and that A and 
A*-^ are subsets of [n]. It will be useful to work with the bipartisation M, A of H, A^ and 
A*-^ which we define as follows. Let m = 2n and let M be the m x m matrix defined 
by the following equations for i,j S [n]: Mij = 0, Mi^n+j = Hij, Mn+ij = Hj^i, and 
Mn+i^n+j = 0. The matrix M can be broken into four "tiles" as follows. 

Let A = A^ U {n + J I J G A*-^}. Note that the matrix Im;A can be decomposed naturally 
in terms of the tiles /^-aa and /„.ac • 

T - ( ^n;AR 

-'m;A — I p. J 

We identify a set of conditions on H, A^ and A*-^ that determine whether or not the 
problem EVAL(M, /m,/m;A) can be computed in polynomial time. We will see how this 
implies Theorem 11.21 

The Group Condition. For an n x n matrix H and a row index / G [n], let 

G{H, I) := {F,,, o Hi^^ I i G [n]} U {-Hi^^ o Hi^^ \ i G [n]} . 

The group condition for H is: 

(GC) For ah / G [n], both G{H,l) = G{H, 1) and G{H^,l) = G{H^, 1). 

The group condition gets its name from the fact that the condition implies that G{H, I) is 
an Abelian group (see Lemma lT.ip . As all elements of this group have order 2, the group 
condition gives us some information about the order of such matrices, as the following 
lemma (which we prove later in Section [7]) shows: 

Lemma 2.1. Let H be an n x n Hadamard matrix. If H satisfies (GC) then n = 2^ for 
some integer k. 



The Representability Conditions. We describe Hadamard matrices H satisfying (GC) 
by F2-polynomials. By Lemma [2. II these matrices have order n = 2^^ . We extend our notion 
of "symmetric representation": Let /) : Fg — > [n] and p ■.¥2 ^ [n] be index mappings 
(i.e. bijective mappings) and X = (Xi,...,Xfc) and Y = {Yi^ . . . ,Yk). A polynomial 
h{X, Y) over F2 represents H with respect to p and p if for all x, y G Fg it holds that 

/i(x,y) = 1 ^^ i^pfl(x),pC(y) = -1- 

So a symmetric representation is just a representation with p = p . We say that the 
set A is linear with respect to p if there is a linear subvectorspace L C F2 a such that 
p^[L^) = A^. Note that, if A^ is linear, then |A^| = 2' for some / < k. We may therefore 
define a coordinatisation of A (with respect to p ) as a linear map <f) : Fg ^ F2 such that 
(/)^(F2) = L^, that is A^ is just the image of the concatenated mapping p^o^-^. We define 
the notion of linearity of A*-^ with respect to p'-' and the coordinatisation of A*-^ with respect 
to p similarly. For a permutation vr E 5^ we use the shorthand X.„ -Y := 0j=x -^iTii) ' ^i- 
The following conditions stipulate the representability (R) of H by F2-polynomials, 
the linearity (L) of the sets A and A , and the appropriate degree restrictions on the 
associated polynomials (D). 

(R) There are index mappings p : ¥2 -^ [n] and p : F2 -^ [n] and a permutation 
vr G Sfe such that (w.r.t. p and p ) the matrix H is represented by a polynomial 
of the form 

h{X,Y) = X^-Yeg'^{X)eg^{Y). (2.1) 

Moreover, if A is non-empty, then p (0) G A . Similarly, if A is non-empty, 

then p'-^{0) G A*-". Finally, if H is symmetric and A^ = A*^, then g^ = g'^ and 
pR ^ pC _ 

(L) A and A are linear with respect to p and p respectively. 

(D) Either A^ is empty or there is a coordinatisation <f)^ of A^ w.r.t p^ such that the 
polynomial g ocf) has degree at most 2. Similarly, either A is empty or there is 
a coordinatisation (j)'~" of A*-^ w.r.t p^ such that the polynomial g^ o cj)^ has degree 
at most 2. Finally, if H is symmetric and A = A is nonempty then (/> = (f) . 

Actually, it turns out that condition (D) is invariant under the choice of the coordinatisa- 
tions (j)^,(j)'-". However, the conditions are not invariant under the choice of the represen- 
tation p ,p , and this is a major source of technical problems. 

Before we can apply the conditions (R), (L) and (D) we deal with one technical issue. 
Let if be an n X n Hadamard matrix and let A , A C [n] be subsets of indices. Let M, A 
be the bipartisation of H, A^ and A*-" . We say that H is positive for A^ and A*-^ if there is 
an entry Htj = +1 such that (1) i e A^ or A^ = 0, (2) j e A'^ or A*^ = 0, and (3) If H is 
symmetric and A = A then i = j. Otherwise, note that —H is positive for A and A . 
Since Zm,i^,i^.JG) = {-1)\''^^^\z_m,i^,i„,.JG), the problems EVAL(M,/„, /^;a) and 
EVAL(— M, /m,/m;A) have equivalent complexity, so we lose no generality by restricting 
attention to the positive case, which is helpful for a technical reason. 

We can now state the theorem which is proved in this section. 



Theorem 2.2. Let H be an n x n Hadamard matrix and let A , A C [n] be subsets 
of indices. Let M, A be the bipartisation of H, A^ and K^ and let m = 2n. If H 
is positive for A and A then EVAL(M, I^, /m;A) is polynomial-time computable if, 
and only if, H A^ and A*^ satisfy the group condition (GC) and conditions (R), (L), 
and (D). Otherwise EVAL(M, I^, /m;A) is ^P-hard. If H is not positive for A and 
A^ then EVAL^M, Im, Im;A) is polynomial-time computable if, and only if, —H A^ and 
A satisfy the group condition (GC) and conditions (R), (L), and (D). Otherwise 
EVAL(M, Im, /m;A) is ^F-hard. There is a polynomial-time algorithm that takes input 
H, A and A and decides whether EVAL(M, /m,,/m;A) is polynomial-time computable or 
ifP-hard. 

The theorem is proved using a sequence of lemmas. Proof sketches of these lemmas will 
be given in this section and full proofs will be given later in Section [71 

Lemma 2.3 (Group Condition Lemma). Let H be annxn Hadamard matrix and let 
A^, A*^ C [n] be subsets of indices. Let M, A be the bipartisation of H , A^ and K'^ and 
let m = 2n. If H does not satisfy (GC) then 'EYAh{M , Im, Im;A) is i^P-hard. There is a 
polynomial-time algorithm that takes determines whether H satisfies (GC). 

Proof sketch. For any integer p and a symmetric non-negative matrix C^^', which depends 
upon H, the proof uses gadgetry to transform an input to EVAL(C^^^) into an input to 
EVAL(M, Im, /m;A)- The fact that H does not satisfy (GC) is used to show that, as long 
as p is sufficiently large with respect to M, then C'-^' has a block of rank greater than 
one. By a result of Bulatov and Grohe, EVAL(CW) is #P-hard, so EYAL{M,Irn,Im;A) is 
#P-hard. 

Lemma 2.4 (Polynomial Representation Lemma). Let H be an n x n Hadamard 
matrix and A^, A^ C [n] subsets of indices. Suppose that H satisfies (GC) and that H is 
positive for A^ and A^ . Then the Representability Condition (R) is satisfied. There is a 
polynomial-time algorithm that computes the representation. 

Proof sketch. The representation is constructed inductively. First, permutations are used 
to transform H into a normalised matrix H, that is, a Hadamard matrix H whose first row 
and column consist entirely of +ls, which still satisfies (GC). We then show that there is 
a permutation of H which can be expressed as the tensor product of a simple Hadamard 
matrix (either H2 or H/^ and a smaller normalised symmetric Hadamard matrix H' . By 
induction, we construct a representation for H' and use this to construct a representation 
for the normalised matrix H of the form X-,^ ■ Y for a permutation tt £ Sk- We use this to 
construct a representation for H. 

Lemma 2.5 (Linearity Lemma). Let H be annxn Hadamard matrix and A^,A^ C [n] 
subsets of indices. Let M, A be the bipartisation of H, A and A and let m = 2n. Suppose 
that (GC) and (R) are satisfied. Then the problem EVAL(M, I^, Im;A) is i^P-hard unless 
the Linearity condition (L) holds. There is a polynomial-time algorithm that determines 
whether (L) holds. 



Proof sketch. For a symmetric non-negative matrix C, which depends upon H, the proof 
uses gadgetry to transform an input to EVAL(C, Im-,Im;A) to an input of EVAL(M, Im, Im-^A.)- 
By (R), there are bijective index mappings /? : Fg — > [n] and p : Fg ^ [n] and a permu- 
tation -K & Sk such that (w.r.t. p^ and p*^) the matrix H is represented by a polynomial 
of the appropriate form. Let r be the inverse of p and r be the inverse of p . Let 
L^ = T^{hC) and L^ = r^(A^). We show that either EVAL(C, /,„,/,„; a) is #P-hard or 
(L) is satisfied. In particular, the assumption that EVAL(C, /m,/m;A)- is not #P-hard 
means that its blocks all have rank 1 by the result of Bulatov and Grohe. We use this 
fact to show that L^ is a linear subspace of A^ and that L'-' is a linear subspace of L^ . 
To show that L^ is a linear space of A^, we use L^ to construct an appropriate linear 
subspace and compare Fourier coefficients to see that it is in fact L^ itself. 

Lemma 2.6 (Degree Lemma). Let H be an n x n Hadamard matrix and A^, A*^ C 
[n] subsets of indices. Let M, A be the bipartisation of H, A^ and A*^ and let m = 
2n. Suppose that (GC),(R) and (L) are satisfied. Then EVAL(M, /mi-fm;A) ^s #P- 
hard unless the Degree Condition (D) holds. There is a polynomial-time algorithm that 
determines whether (D) holds. 

Proof sketch. For any (even) integer p and a symmetric non-negative matrix C'-'^' , which 
depends upon H, the proof uses gadgetry to transform an input to EVAL(C^^^) into an 
input to EVAL(M, Jm, /m;A)- Using the representation of H, a coordinatisation (p^ with 
respect to A , and a coordinatisation (j) with respect to A , some of the entries Cj^l of 

the matrix C'-^' may be expressed as sums, over elements in Fg, for some i, of appropri- 
ate powers of —1. We study properties of polynomials g{Xi, . . . ,Xk) G F2[Xi, . . . ,Xk], 
discovering that the number of roots of a certain polynomial ga,f3,'yiXi, ■ ■ ■ , X].), which is 
derived from g{Xi, . . . , Xk), depends upon the degree of g. From this we can show that if 
(D) does not hold then there is an even p such that EVAL(C[p1) is ^^P-hard. 

Proof (Proof of Theorem \2.2\) . By the equivalence of the problems EVAL(M, IrmIm;A) and 
EVAL(— M, Lm, Im;A) we can assume that H is positive for A and A . The hardness part 
follows directly from the Lemmas above. We shall give the proof for the tractability part. 
Given H, A and A satisfying (GC), (R), (L) and (D), we shall show how to compute 
ZujmJmAi^) ^^^ ^^^ input graph G in polynomial time. 

Note first that ZMjrnJmAi^) ~ ^ unless G is bipartite. If G has connected components 
Gi, . . . Gc, then 

c 
Zm,I^J^,a{G) = ||-^M,/„,7„;A(Gi). 

«=1 
Therefore, it suffices to give the proof for connected bipartite graphs. Let G = {V, E) be 
such a graph with vertex bipartition U (JW = V . Let 1^ ^ ^ be the set of odd-degree 
vertices in G and let Uo = W f\Vo and Wo = W f\Vohe the corresponding subsets of U 
and W. Let Ue = U \Uo imdWe = W \ Wo- We have 

i:V^[m] {u,w}&E vGVo ^iV^lm] {u,w}eE 

?(Vo)CA 



As G is bipartite and connected this sum splits into ZM,i,n,i^./^iG) = Z^ + Z^ for values 



z- 



E E n ^. 

C-U^[n] (-.W^ln] {u,w}eE 



5(«),CM and Z^= ^ ^ JJ 

i-.U^ln] C-W^ln] {u,w}eE 



i?, 



C(«'),5(«) 



We will show how to compute Z'^ . The computation of the value Z^ is similar. 

Fix configurations $, : U ^ [n] and (^ : VF ^ [n] and let p^, p*^ be the index mappings 
and h the F2-polynomial representing H as given in condition (R). Let r be the inverse 
of p^ and let r'^ be the inverse of p'^. Let L^ = t^{A^) and L*-" = r'-"(A'^). Then ^ and 
( induce a configuration ^ : y ^ Fg defined by 



<i{v) :-- 






which implies, for all u (i U,w (i W that /i(^(u),?(ttj)) = 1 iff H^fy_\^ty^\ = —1. Let (j)^ and 
(/) be coordinatisations of A and A w.r.t. p and p satisfying (L) and (D). We can 
simplify 



2- = E E n ( 

(,:U^[n\ C-W^ln] {u,w}gE 

aUo)cAR C{Wo)cac ueu 

= E 



-1 



,A(r«(5(«)),rC(C(«,))) 



-1 



){n,»}ei5:«e!7 '*('?(«)'?("')) 






Define, for a G F2, sets 



? : y ^ F^ I <j(?7o) C L^, <j(Ty„) C L^, h{^{u),^{w)) = a 



{u,w}eE 



(2.2) 



Then Z^ = sq — si. Therefore, it remains to show how to compute the values Sa- Define, 
for each w G F, a tuple X"" = {Xf, . . . , XJl) and let he be the F2-polynomial 



ho- h{X-,X- 



{u,w}£E 



{X^\ • X- © g^{X^) © g^{X^). (2.3) 



{u,w}eE 
ueu 



Here the second equality follows from the definition of the polynomial h given in condition 
(R) and the fact that the terms g^{X^) and g'--' {X"^) in the definition of h appear exactly 
deg(u) and deg(tt;) many times in ho- Therefore, these terms cancel for all even degree 
vertices. 

Let var(/iG') denote the set of variables in he and for mappings x '■ var(/iG) -^ F2 we 
use the expression x(X^) := {x{Xi), . . . ^x{^k)) as a shorthand and define the F2-sum 
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hcix) ■= ®{u,w}&E:ueuKx(.X''),x(.X'")). We find that Sa can be expressed by 

n \ ml xiX"") £ L^ for all -u G C/o, ,, . .1 ,^ ,. 

^ ■ ^^^("«) - ^' I xiXn G LO for all w G T^J, ^^^^ = "^j ^^''^ 



By equation (12. 4p we are interested only in those assignments x of the variables of 
ha which satisfy x(^") e ^^ and x(^"') e L'^ for all u G [/„ and w G Wo- With 
|A^| = 2^ and jA*^! = 2^ for some appropriate £^,i'-^, we introduce variable vectors 
y" = (y«, . . . , Y^« ) and Z"" = {Zf, ..., Z^c) for all u e Uo and w e Wo- li u e Uo or 
w G Wo then we can express the term {X'^)„ ■ X^ in ho in terms of these new variables. 
In particular, let 

hh = i<P''{Y^)).-^^iZ^)(B {Xn.-X^ 

{u,w}eE {u,w}eE 

ueUo,w€Wo ueUe,w£We 

{u,w}<^E {u,w}eE 

ueUe,WGWo ueUo,WGWe 

Let 

h'c = /^G © 5^(<^''(^")) © g'^i^P'^iZ-)) (2.5) 

We therefore have 

sa = \{x- var(/i'G) ^ F2 | h'^ix) = a)}\. (2.6) 

By condition (D), the polynomials g o cf)"- and g o (jf' are of degree at most 2 and 
therefore H'q is a polynomial of degree at most 2. Furthermore, we have expressed Sa as the 
number of solutions to a polynomial equation over F2. Therefore, the proof now follows 
by the following well-known fact. 

Fact 2.7. The number of solutions to polynomial equations of degree at most 2 over F2 
can he computed in polynomial time. 

This is a direct consequence of Theorems 6.30 and 6.32 in [19j (see also [13]). □ 

3 The General Case 

In this section we will prove Theorem 11.11 Before we can give the proof some further 
results have to be derived, which then enable us to extend Theorems 11.21 and 12.21 It 
will be convenient to focus on connected components. This is expressed by the following 
lemma, which will be proved later in Section [8l 

Lemma 3.1. Let A he a symmetric real-valued matrix with components Ai, . . . , A^. Then 
the following holds 

(1) //EVAL(^i) is iJ^V -hard for some i G [c] then EVAL(A) is #P-hard. 
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(2) If EVAL(^i) is PTIME computable for all i G [c] then EVAL(yl) is PTIME com- 
putable. 

Recall that for each connected symmetric matrix A there is a block B such that either 

A = B or, up to permutation of the rows and columns, A = ( j^ J . We call B 

the block underlying A. For such connected A we furthermore see that the evaluation 
problem is either ^P-hard or we can reduce it to the evaluation problem on bipartisations 
of Hadamard matrices. This is expressed in the following lemma, which will be proved 
later in Section [51 

Lemma 3.2. Suppose that A is a symmetric connected matrix. Then either EVAL(yl) is 
f^V-hard or the following holds. 

(1) If A is not bipartite there is a symmetric rxr Hadamard matrix H and a set A C [r] 
such that 

EVAL(^) = EYAL{H,Ir,Ir;AR)- 

(2) If A is bipartite then there is an r x r Hadamard matrix H, sets A^, A*-^ C [r] and a 
bipartisation M, A of H, A and A such that 

EVAL(^) = EVAL(M, hr, hr.K)- 

Furthermore it can be decided in time polynomial in the size of A which of the three 
alternatives (^V -hardness, (1), or (2)) holds. 

We are now able to prove the main Theorem. 

Proof (Proof of Theorem \l.l\) . Given a symmetric matrix A E M™^™. By Lemma |3. II we 
may assume that the matrix A is connected. By Lemma 13.21 Theorem 12.21 and Corollary 
17.141 the problem EVAL(yl) is either polynomial time computable or ^^P-hard. The ex- 
istence of a polynomial time algorithm for deciding which of the two possibilities holds, 
given a matrix A, follows directly by these results. □ 

4 Outline of the paper 

The rest of the paper is organised as follows. Section [5] describes some generalisations 
of the partition function evaluation problem, including related work and open problems. 
The remaining sections contain the proofs of the lemmas which have been stated without 
proof. In particular. Section [6] develops some technical tools which we will use. Section [7] 
proves the lemmas that are stated in Section [2l Finally, Section [8] proves the lemmas that 
are stated in Section [3l 
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5 Related Work 

There are several natural directions in which the work in this paper could be extended. 
First, the matrix A could be extended to include algebraic complex numbers and not merely 
algebraic real numbers. This extension has been carried out, subsequent to this paper, in 
an impressive 111-page paper by Cai, Chen and Lu [8j. The work could also be extended by 
allowing the matrix A to be asymmetric. A recent breakthrough by Bulatov [7j establishes 
the existence of a dichotomy theorem for counting satisfying assignments in constraint 
satisfaction. This implies that a dichotomy exists for the case in which A is a 0-1 matrix 
(which is not necessarily symmetric) — in this case computing the partition function 
corresponds to counting homomorphisms to a directed graph, in particular to the directed 
graph with adjacency matrix A. Bulatov's dichotomy is not known to be effective in the 
sense that it is not known to be decidable given a matrix A whether it is tractable or 
not. An effective dichotomy was given by Dyer, Goldberg and Paterson [11] for the special 
case in which the directed graph with adjacency matrix A is acyclic but no effective 
dichotomy is currently known for the case of general 0-1 matrices A. A generalisation of 
Bulatov and Grohe's dichotomy for symmetric non-negative matrices to symmetric non- 
negative functions of arbitrary arity was given recently by Dyer, Goldberg, and Jerrum 
[lOj . However, nothing is known about how how to handle functions of larger arity in the 
presence of mixed signs. 

6 Technical Tools 

6.1 Stretchings and Thickenings 

We introduce some fundamental relations which will be used in most of our reductions. 
Let G = (y, E) be a graph. The s-stretch of G is the graph SgG obtained from G by 
replacing each edge by a path on s edges. The t-thickening of G is the graph TtG obtained 
from G by replacing each edge by t parallel edges. Let A^^' denote the matrix obtained 
from A by taking each of its entries to the power of t. 

Lemma 6.1 ([9j). For a symmetric matrix A G M^^™ and a diagonal m x m matrix D 
we have, for a// s, t G N 

EYAL{A{DAY-^,D) < EVAL(^,L>) and EYAL{A^^\ D) < EYAL{A,D) 

These reducibilities hold as 

^A{DAy'^,DiG) = Za,d{SsG) and Z^(t) £,(G) = Za,d(TjG). 

6.1.1 Twin Reduction 

We need some extensions of Lemma 3.5 in [9]. For a symmetric m x m matrix A we say 
that two rows Ai^^ and Aj^^ are twins iff Ai^^, = Aj^^,. This induces an equivalence relation 
on the rows (and by symmetry on the columns) of A. Let /i, . . . /„ be a partition of the 
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row indices of A according to this relation. The twin- resolvent of A is the matrix defined, 
for ah i,j G [n], by 

T{A)ij := A^^i, for some fi £ li,^ £ Ij. 

The definition of the classes /j implies that A^^i, = A^/y for all fi, fi' G /j and z^, u' G Jj 
and therefore the matrix T(A) is well-defined. 

The above definition furthermore give rise to a mapping r : [m] -^ [n] defined by 
/i G ^T(ri) that is T maps /j, G [m] to the class Ij it is contained in. Therefore, we have 
T{A)^u\^rj\ = Aij for all i,j G [m]. We call r the twin-resolution mapping of ^. 

Lemma 6.2 (Twin Reduction Lemma). Let A be a symmetric rax m matrix and D 
a diagonal mx m matrix of vertex weights. Let /i, . . . , /„ be a partition of the row indices 
of A according to the twin-relation. Then 

Za,d{G) = Zr(A),A(^) f"'^ "^^ graphs G 
where A is a diagonal n x n matrix defined by Aj^j = X^^g/. -Dj/,;/ for all i G [n]. 
Proof. Let r be the twin-resolution mapping of A. Then 

za,d{g) = X] n ^?(«).5w n ^?(^)-c(^) 

^■.V-^[m] {u,v}eE veV 

= Z] n '^(^)ro?(«),ro5{«) n ^C(^):C(^) 

where the second equality follows from the definition of r. As for all ^ : y ^ [m] we 
have T o ^ : V ^ [n], we can partition the S, into classes according to their images under 
concatenation with r and obtain: 

Za,d(.G) = 22 z2 li '^("^)V'(«),V'M II -^5(«),5{«) 

V':V-»[n] g:y^[m] {u,v}eE v€V 

I \ 



ip:V^[n] {u,v}gE 



\ T0S,=1p I 



Fix some il^ : V ^ [n]. For ^ : V ^ [m] we have r o ,f = ^ if and only if ^ ({^}) = C i^i) 
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for all i G [n]. Define Vi := ip ^{{i}) for all i G [n] which yields a partition of V. Thus 
Y. n ^^yUiy) = Z] n ^^yUiy) 

Toi=iP \^ ie[n]: ^(V^)Ch 

n 

= n iz n ^ii^uH 

i='i-^i:Vi^IiVeVi 
n 

= nnEo- 



Hence 



^■.V^[n]{u,v}eE veV 



6.2 Basic Tractability and #P-hardness 

The following Lemma is a straightforward extension of Theorem 6 in [6] . 

Lemma 6.3. Let A E M^^™ be a symmetric m,atrix and D a diagonal m x m matrix. If 
each component of A either has row rank 1 or is bipartite and has rank 2 then EVAL(j4, D) 
is polynomial time computable. 

Proof. Let G = {V, E) be a given graph with components Gi, . . . , Gc and let Ai, . . . .,Ai 
be the components of A and Di , . . . ,Di the submatrices of D corresponding to these 
components. Then 

c I 

ZA,D{G) = Y{Y.^^.MGi)- 
i=l j=l 

Therefore the proof follows straightforwardly from the special case of connected G and A. 
Assume therefore that both G and A are connected. 

We will prove the following claim, which holds for directed graphs. 

Claim 1. Let B"^^"^ be a (not necessarily symmetric) matrix of row rank 1 and D' a 
diagonal matrix. Then for every directed graph G the value 

^■.V^[m] {u,v)£E veV 

can be computed in polynomial time. 
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Proof. Let G = {V, E) be a directed graph and for every vertex v £ V denote by outdeg(f ) 
and indeg(i') the number of outgoing and incoming edges incident with v. There are vectors 
a,b £ W^ such that B = ab . Then, for every configuration ^ -.V ^ [m], 

no TT J, TT outdog(D),indeg(ti) 



{u,v)&E {u,v)£E v£V 



and therefore 



^■.V^[m] {u,v)eE veV 

ETT outdeg{i;),indeg(i)) ^/ 
11 "?(«) '^av) ^iivuiv) 

C.V^[m]v(^V 

m 
_ TT- V^^outdeg(D)^indeg(t;)^/ 

v&V i=l 

And the terms in the last hne can be evaluated in polynomial time. This completes the 
proof of the claim. 

With this claim we are now able to prove the Lemma. Recall that A is connected and 
symmetric. If A is non-bipartite then A has rank 1. For a given connected graph G 
let G' be a directed graph obtained from G by orienting its edges arbitrarily. We have 
Za,d{G) = Z*j^j^{G') and the value Z\j-){G') can be computed by Claim [TJ 

Otherwise, if A is bipartite then we have (up to permutation of the rows/columns of A) 

A-( ° ^ 
for a block B of rank 1. Let A' be the matrix 

^ - Vo 

which has rank 1 because B has. Note furthermore, that Za,d{G) = unless G is bipartite. 
Assume therefore that G = {U,W,E) is a bipartite graph and let the graphs Gjjy/, Gwu 
be obtained from G by directing all edges from U to W {W to U, resp.). Then 

Aa,d{G) = Za',d{Guw) + Za',d{Gwu) 
and the terms of the right hand side are polynomial time computable by Claim [H □ 

The following ^^P-hardness result will be the basis of all our proofs of intractability. 

Lemma 6.4. Given a symmetric matrix A of order n and diagonal n x n matrices D, O 
such that D is a non-singular matrix of non-negative integers. If ahs{A) contains a block 
of row rank at least 2 then EVAL(yl, D, O) is #P-hard. 
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Proof. Observe that by 2-thickenmg we have EYAL{A^^\D) < EYAL{A, D,0). We can 
form a matrix A' from A^"^' by introducing twins according to D that is, doing the inverse 
operation of Lemma I6.2[ More precisely, let rii := Di^i for all i G [n] and define m := 
n ■ {Y17=i ^«)- '^° define the m x m matrix A' we consider its row and column indices as 
pairs and define 

^Ui),{X,j) '■= ^S for all K,\G[n\,ie n^, j G nx. (6.1) 

By the definition of A' we see that Application of the Twin Reduction Lemma 16.21 to A' 
yields 

Za'{G) = Z^(2) £)(G) for every graph G. 

and thus EVAL(A') = EVAL(^(^), D). By equation ([6?T]) the matrix A' contains a block 
of row rank at least 2 iff A^"^' does which in turn is the case iff abs(^) contains such a 
block. The proof now follows from the result of Bulatov and Grohe |6j. □ 

6.3 Interpolation Lemma 

In the next chapters we will make extensive use of the following lemma which is an analogue 
of the interpolation technique as used for example in [9j . 

Lemma 6.5. Let xi, . . . ,Xn G lR>o be pairwise distinct and let V and N two finite multi- 
sets of real numbers with \V\ = \M\ = n. Then the following are equivalent 

(1) V=M 

(2) there is an ordering of the elements in V and J\f such that for arbitrarily large p, we 
have 

Proof. The forward direction is trivial. Hence, assume that (2) holds but not (1). With 
the given ordering of V and M we have V = {ai, . . . , a„} and M = bi, . . . ,bn. We may 
assume that there is no z S [n] such that Oj = bi because otherwise, we might delete this 
pair from V and M. Hence, let k £ [n] be such that x^ = maxjg[„] Xj. Assume w.l.o.g. 
that ak > bk then, for a constant c 7^ 

= ^ x^ai - Y ^^^i = ^ki^-k -bk) + Y ^^^"'^ ~ ^*) 

ai&V bi£M ie[n]\{fc} 

ie[n]\{fc} ^ '^^ 



\ p 

Xi 



By hmp_>oo Y^i(p\n]\^k\ \ — ) i'^i ~ ^i) = 0' ^^^^ yields a contradiction. □ 
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Lemma 6.6. Let xi, . . . ,Xn £ M>o be pairwise distinct and /ei ai, . . . , a„ G M and bi, . . . ,bn G 
M. There is a po €f^ such that for all p > po, the equation 

n n 

1=1 1=1 

holds if, and only if, ai = bi for all i £ [n]. 

Proof. Note first that backward direction is trivial. It remains therefore to prove the 
fohowing. For each / C [n] there is a p/ G N such that for all p > pi, if 

^x^,ia,-bi) = (6.2) 

then Oj = bi for all i £ I. We will give the proof by induction on the cardinality of /. For 
empty / there is nothing to be shown. Assume therefore that / 7^ let A: G / be such that 
Xk = maxjg/ Xi and define I' = I\{k}. 



Claim 1. There is a p^ G N such that for all p > p^, ii equation (j6.2p is satisfied then 
a-k = bk- 

Proof. Assume for contradiction that a^ 7^ bk but equation (j6.2p holds for all p G N. This 
implies 

= {ak-bk) + Y.(-y{a,-h). (6.3) 

iei' ^^^/ 

As i £ I' with Oj = bi do not contribute to the above sum we may further assume that 
ai 7^ hi for all i £ I' . If /' = we already have a contradiction. If otherwise /' 7^ 0, let k' 
be such that xy = maxjg// Xj. We find that 



Xi 



SiSX 



— y \ai-bi\. 



< 

\ XL. I 

p 



In particular equation (j6.3p does not hold if I ] ^jg// |aj — 6j| < \ak — bk\ which, as 

\ ^k J 
Xk > ^k' is the case for all 

log \ak - bk\ - log Y^idV l«i - ^i\ 



p> 



{\ogXk' -logxk) 
in contradiction to our assumption. H 

By the induction hypothesis there is a p/' such that for all p > pp 

Y^x^{a,-bi) = 

implies aj = 6j for all i £ I'. Let pk be defined as in Claim 1 then the proof follows with 
pj = max{pk,pj,}. n 
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7 The Proofs for Section [2] 

7.1 Notation and Preliminaries 

For X = {xi, . . . , Xn),y = (yi, . . . , yn) £ I^"i by (x, y) we denote the inner product Y17=i ^iVi 
of X and y. It may be a source of confusion that we work over two different fields, M and 
F2. Addition in F2 is denoted by ©, and for a = (ai, . . . , a„),/3 = (/3i, . . . ,/?„) G Fg, 
a • /? is the dot product 0j=iai/3i. Similarly, for vr G 5*^, ajr • /? denotes 0j=i a7r(^)/9^• 
a © /? denotes the element (ai © /3i, • • • , a^ © /3fc) in Fg. Similarly, for vr G S'fc, Ott © /9 
denotes the element (a7r(i) © A, • ' ' > cn-K(k) © /^fc)- Similar notation applies to variables, so 
if X = (Xi, . . . , Xfe) and y = (Fi, . . . , Yk) then X^ • F denotes 0^^^ X^(^i)Yi. For / C [fc], 
let Let X \ / be the tuple containing, in order, all variables in {Xi, . . . ,Xfc} other than 
those with indices in /. For example, X \ {2, 3} denotes the tuple (Xi, X4, . . . , X^). 

7.2 Tiie Group Condition 

Lemma 7.1. Let H he an n x n Hadamard matrix. If H satisfies (GC) then G{H, 1) 
forms an Abelian group under the Hadamard product. 

Proof. Commutativity and associativity follow from the definition of the Hadamard prod- 
uct. To show closure, we consider two elements in G{H, 1) and show that their Hadamard 
product is also in G{H,1). First, consider //j^, o i^i,^ and Hj,^ o ffi,^. Their Hadamard 
product is Hi.^ o Hi^:^ o Hj:^ o Hi^^ = Hi^^ o Hj^^ which is in G{H,j) by the definition of 
G{H,j) and therefore in Hi^^ by (GC). Similarly, we fined that the product of —Hi^^^oHi^^, 
and Hj^, o Hi^^, is in G{H, 1) and also the product of —Hi^^, o Hi^, and —Hj^^ o Hi^^, is in 
G{H, 1). From closure, it follows that the product of Hi^^ o Hi^^ and itself is in G{H, 1) 
and this row (the all ones row) is the identity element in the group. □ 

Proof (Proof of Lemma \2.1\) . By Lemma |7. 11 G{H, 1) forms an Abelian group under the 
Hadamard product. All elements of this group have order 2, and thus it follows from 
elementary algebra that the order of the group is a power of 2. □ 

Proof (Proof of Lem,ma \2.3[ the Group Condition lemma). It is clear from the definition 
of the Group Condition that there is a polynomial-time algorithm that determines whether 
H satisfies (GC). We focus on the #P-hardness result. Let EVALeven(yl) denote the 
problem of computing Zyi(G) for an input graph G in which every vertex of G has even 
degree. 

Let H, n, M, A and m be defined as in the statement of the lemma. Let p be an 
even number. We will show how to transform any graph G into a graph Gp with all 
even-degree vertices so that Z(j[p]{G) = Zm{Gp) for a matrix C'^' which we will define 
below. The definition of Cl^l depends upon M but not upon G. Thus, we will have 
EVAL(CW) < EVALeven(M) < EVAL(M,/™,/m;A). 

To finish the proof, we will show that, as long as p is sufficiently large with respect to M, 
then EVAL(CW) is #P-hard. 
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Figure 1: The gadget for p = 4 
We start by giving the transformation from G = (y, E) into Gp = (Ip, Ep): 

Vp := V U {Ve,Ve'^,Ve,l... ,Ve,p\e £ E} 

Ep := {{u,Ve,i} ,...,{u,Ve,p} \e = {u,v} £ E} 
U { {v, Ve,l} ,■■■ ,{v, Ve,p} \e = {u,v} G E} 
U { {Ve,l,Ve} ,..., {Ve,p, fe} | 6 G E} 
U {{Ve,l,Ve''} , . . . ,{Ve,p,Ve'^} \ e £ E} 

EssentiaUy, every edge e = {n, v} in G is replaced by a distinct gadget. Figure [U illustrates 
this gadget for p = 4. Since p is even, it is clear that all vertices of Gp have even degree. 

Let us now construct the matrix G^'^'. Let F denote the graph with vertices u and v 
and a single edge between them. Clearly C^ ■ is equal to the contribution to Zm^^^) 
corresponding to those configurations S, with ^{u) = i and ^{v) = j. Thus, 



m m 



^S = E E E M.,cM,,,Ma,cMb,c , (7.1) 



a=l 6=1 \c=l / 

where a denotes the choice of spin for Ve and b denotes the choice of spin for Ve^ and c 
denotes the choice of spin for a vertex Ve/- 

To finish the proof we must show that, as long as p is sufficiently large with respect 
to M, then EVAL^G^^) is #P-hard. From the definition of M, we see that, for i £ [n], 
j £ {n + 1, . . . , 2n}, we have C^ ■ = G^- = 0. Also, for all i,j £ [n], we have the following. 



eg = E E ^^i'* ° ^J.*' ^-.* ° Hb,*f, and 

a=l b=l 



Cn+i,n+j - E E ^^*'^ ° H*,j^H*,a ° H^ 



a=l b=l 
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Now, for all i,j € [n] and x G {0, . . . , n} let sf^ be the number of pairs (a, b) such that 

|(i/j^* o ifj^*, -ffa,* o -f^6,*)l = 2; and similarly let s|^_[_j ^^ • be the number of pairs (a, b) such 
that |(^*,i o H^,j,Ht:^a o -ff*,6)| = X. Then for all z,j G [n] we have 

?i n 

rylp] _ V^ „N p J ^H _ V^ „H ^p (J ON 

«J ~ Z^ hj '-■n+i,n+j — Z^^n+i,n+j^ ■ U'^J 

x=0 x=0 

The pair (a, 6) = (i,j) contributes one towards s^"- and one towards s^^-^_^- so, for all 
^) J £ i'^]) we have CJ'^J > and Cjf_|_^ ^^ • > (remember that p is even). 

Since H is Hadamard, s^" = n for every i S [n] and, for every x G {1, . . . , n — 1}, s^^ = 
so dfj = riP^^. Also, since H is Hadamard, HH^ = nl, so H^ /n is the right inverse, 
hence also the left inverse, of H, so {l/n)H^ H = /, so H^ is also Hadamard. It follows 
that sj![i „+, = n and, for every x G {1, . . . , n - 1}, s^^li^^^^ = so c|fj, „+i = nP+^ 

We will prove that EVAL(Cl^]) is #P-hard for some sufficiently large even p. We will 
assume for contradiction that, for every even p, EVAL(C'^') is not ^^P-hard. Equation (j7.ip 
indicates that C^^^ is symmetric, so by Lemma 16.41 (due to Bulatov and Grohe), for every 
even p, both blocks of C^^^ have rank 1. This means that every principal 2x2 submatrix 
in the blocks has a zero determinant. So, for i,j G [n], we have (C|^ )^ — {Cfh'^ = and 

i^n+i,n+i) ~ (C'n+i,n+j) = 0, SO 

Since equations (j7.2p and (|7.3|) hold for all even p and all i,j G [n]. Lemma 16.51 allows us 
to deduce that, for all i,j G [n] and x G {0, . . . ,n}, sfj = sfj and s„|j„+j = s|^|j^„+j. 
Thus, for all i,j G [n], 

Jl] _ . . . _ „["-!] _ „[1] - ... - J"-^] - n arid ^["1 - ^["1 - n C? 4^ 

From the statement of the lemma, we assume that H does not satisfy (GC). There are 
two similar cases. 

Case 1: Suppose there are i,j G [n] such that G{H,i) ^ G{H,j). Fix such a pair i,j. 
Fix a G [n] such that ii^a,*o-ffj,* is not in G{H,j), Now consider any 6 G [n]. If it were the 
case that \{Ha^<toHi^^, -ff6,*o^j,*)| = n then we would know that either Ha^vHi^y = Hb^^Hj^y 
for all V or Ha,vHi^y = —H^^yHj^v for all v. Either of these would imply Ha^*oHi^^ G G{H,j) 
which is not the case. So we conclude that \{Ha^:t o Hi^^, Hh^* o Hj^^)\ < n. 

Furthermore, there is some b G [n] such that |(-ffa,* ° Hi^^^, Hh^^ o Hj^^,)\ / 0. Otherwise, 

would be a set of n + 1 linearly independent vectors, which is impossible. 

But this implies that for some x G [n — 1] we have s\^- ^ contradicting equation (j7.4p . 
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Case 2: Suppose there are i,j G [n] such that G{H ,i) ^ G{H ,j). As in Case 1, we 
can deduce that \{HJ^^ o HJ^,H'^^ o HJ^)\ < n. Furthermore, there is some b S [n] such 
that K-ff^* o Hf^,H^^ o HJ^)\ 7^ 0. But this imphes that for some x G [n — 1] we have 
■s„_l_j_„+j / contradicting equation (|7.4p . □ 

7.3 Polynomial Representation 

For an n X n matrix H and a row index I £ [n], let R{H) := {-ffj,* | i G [n]}. The Extended 
Group Condition for H is: 

(EGC) R{H) is an Abehan group under the Hadamard product. 

The fohowing lemmas are useful preparation for the proof of Lemma [2.4l the Polynomial 
Representation Lemma. We say that a Hadamard matrix is normalised if its first row and 
column consists entirely of +ls. 

Lemma 7.2. Let H he a normalised nx n Hadamard matrix. If G{H, 1) is closed under 
the Hadamard product then R{H) is closed under the Hadamard product. 

Proof. Fix i,j S [n]. Since G{H, 1) is closed under the Hadamard product, and Hi^^^oHi^^ E 
H{G, 1) and Hj^^ o Hi^^ £ H{G, 1), we have Hi^^ o Hj^^, G H{G, 1). Thus, there is a £ G [n] 
such that either Hi^^^oHj^^, = H^,,oHi^^ = H^^^, (using the fact that the first row of H is all 
ones) or i?j_* o fZ^- „, = —Hi^^oHi^^: = —Hi,,. The latter is equivalent to Hi^^oHi^ = —Hj^^,. 
And since Hj^i = 1 (since the first column of H is positive) this implies that one of Hii 
and H^^i is negative, a contradiction. We conclude that Hi^^, o Hj^^ = H^,^. □ 

Corollary 7.3. Let H be a normalised n x n Hadamard matrix. If H satisfies the group 
condition then H satisfies the extended group condition. 

Proof. Suppose that H satisfies the group condition. By Lemma [7T1 G{H, 1) is an Abelian 
group under the Hadamard product. The identity is the all ones row, which is in R(H), and 
every element is its own inverse. Closure of R{H) follows from Lemma 17.21 □ 

Lemma 7.4. Suppose that B is an r x r matrix with entries in {—1,+!} and that C is 
an t X t matrix with entries in {—1, +1}. Suppose that the tensor product H = B <5^ C is 
a Hadamard matrix. Then B and C are Hadamard. If H is symmetric then so are B and 
C. If H and B are normalised and H satisfies (EGC), then B and C satisfy (EGC) and 
C is normalised. 

Proof. Since H is Hadamard, we know that for any such k £ [r] and distinct i and i' in 
[t], the inner product (i:f(fe_i)i+i,*,i^(fc-i)i+j',*) is zero. But this inner product is 

/ , / , -f^(A:-l)t+i,(£-l)t+i-f^(A:-l)t+i',(^-l)t+j ^ Z_^ Z^ Bk/CijBk/Ci'j 

ie[r]je[t] i&[r]jm 

= f\Gi^^, Ci'^=K/, 
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so C is Hadamard. Similarly, for any distinct k,k' G [r] and any i £ [t], 

= {H(^k-i)t+i,*, H(^k'-i)t+i,*) = 2^ 2^ -f^(fc-i)t+j,(£-i)t+j-f^(fc'-i)t+i,(^-i)t+j 

te[r]jG[t] 

te[r]je[t] 

= t{Bk,*,Bk'^*), 

so S is Hadamard. If -ff is symmetric then it is easy to see that B and C are symmetric as 
well. Also, if H and B are normalised, then it is easy to see that C is normalised as well. 

Suppose now that H and B are normalised and H satisfies (EGC). We first show that 
C satisfies (EGC). Then we will finish by showing that B satisfies (EGC). 

To show that R{C) is an Abelian group under the Hadamard product we just need to 
show closure. (Commutativity and Associativity come from the definition of the Hadamard 
product, the identity element is the row of all ones, and every element is its own inverse.) 
Since R{H) is closed under the Hadamard product, we know that, for any distinct i, i' £ [t], 
Hi^^ o Hif^^: G R{H). But the first t elements of this row are HiiHi/i, . . . ,Hi^tHi'^t = 
Bi^iCi^iBi^iCii^i, . . . ,Bi^iCi^tBi^iCi'^i which is equal to Cj_,, o Cj'^*. This shows that Cj^* o 
Cj/^^. G G(C, 1). Now use lemma [7^2] to show that R{C) is closed under the Hadamard 
product. 

Similarly, to show that R{B) is closed under the Hadamard product, note that for any 
distinct k, k' G [r], -ff(fc_i)t+i * o H (j,i _i~^i_^i^^, G R{H). But the elements of this row are 

-f^(fe-l)t+l,(f-l)t+j-f^(A;'-l)t+l,(£-l)t+j) 

for £ G [r], j G [t], and taking those with i = 1 (which occur every r elements along 
the row) we get -Bfc,iCiji?fc/^iCij-. Thus, the sub-row of these elements is the Hadamard 
product of -Bfc,* and -Bfc',*. This shows that i?fc^* o Bf^/^^ G G{B, 1). Now use lemma 7L2\ to 
show that R{B) is closed under the Hadamard product. □ 

Given an n x n matrix H and permutations S and H in S'n, let i/E,n denote the matrix 
with (-ffs,n)ij = -f^s(i),n{i)- 

Lemma 7.5. Let H be a normalised n x n Hadamard matrix with n > 2 that satisfies 
(GC). Then there are permutations E,n in Sn with S(l) = 1 and n(l) = 1 and a 
normalised Hadamard matrix H' satisfying (GC) such that Hy,^yi = H2 © H' . S, H, and 
H' can he constructed in polynomial time. 

Proof. By Lemma l2. II we know n is a power of 2, say n = 2^^. The lemma is trivial for 
/c = since H = H2 and S and H can be taken to be the identity. So suppose A; > 1. 
Let u = 2^ 

Part 1: Choose S' and H' in Sn with S'(l) = 1 and n'(l) = 1 so that {HY,'^Ti')u+i,u+i = 
-1. 
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How to choose S' and IT': H is Hadamard, so some entry Hij = — 1. The indices i and 
j are not 1 because H is normalised. Let S' be the transposition (i, v + 1) and let 11' be 
the transposition (j, z/ + 1). 

Part 2C: Choose vr in 5„ with 7r(l) = 1 and tt{v + 1) = z^ + 1 so that, for £ G [z^], 

(^s',n"),.+i,£ = +1 and (-H'E',n")j.+i,^+£ = -1> (7-5) 

where 11" denotes the composition of first 11' then tt. 

How to choose vr: We construct a sequence of permutations tti,. . . ,tTi, where vri is the 
identity and we let vr = n^. Let H^ denote -ffs'.Tr-n'- For j G {2, . . . z^}, we define vtj as 
follows. If Hl_^^ ^j^- = —1 then ttj = ttj-i. Otherwise, there isanl<£<i^ + l with 

^i+i e ~ "■*-• ^° ^i' ^^ ^^^ composition of first applying 7r'_^ and then transposing u + j 
and £. To see that such an £ exists, note that H is Hadamard, so {Hi^^, H^j^i^^) = 0. But 
Hi,, is positive, so Hyj^i^^ has exactly v ones, i > 1 because 7rj_in'(l) = 1. 

Part 2R: Choose a in 5„ with a{l) = 1 and a{i' + 1) = z^ + 1 so that, for ^ S [v\, 

(^s",n")£,^+i = +1 and (-H'E",n"),.+^,^+i = -1> (7-6) 

where S" denotes the composition of first S' then a. 

How to choose u: This is symmetric to how we chose vr. 

Since a{v + I) = v + 1, we have (-ffE',n"),.+i/ = {H^'J^")a(u+i)/ = (-^S",n")i.+i,^ for 
every £ G [n], so Equations (|7.5p and (j7.6p give 

(i?s",n"),.+i/ = (-f^E",n")^,i.+i = +1 and (i;rs",n")i.+i,,.+<? = (^s",n")i.+£,,.+i = -1. 

(7.7) 
Part 3C: Choose tt' in S'„ with 7r'(l) = 1 and 7r'([z^]) = [z^] so that, for j,l G [z^] 

(-f^E",n)^j = (i^s",n)^^^+j and (i;^s",n)^+^j = -{Hj]",n)^_^(^^^j^p (7.8) 

where H denotes the composition of first H" then vr'. 

How to choose vr': Note that H satisfies (EGC) by Corollary 17.31 hence Hy,"ji" satisfies 
(EGC) (permuting does not change (EGG)). Start with vr'(l) = 1 and vr'(i/ + 1) = z/ + 1. 
Note that for j = 1, we have, by normalisation and Equation (17.61) . 



V£e M,(^s",n")f,^'(j) = (^s",n")£,^'(^+j) and (-f^E",n"),.+f,^'(j) = -(-^s",n")^+^,^'(i.+j): 

(7.9) 
Now for j G {2, . . . , z^} we define vr'(j) and vr'(z^ + j) to satisfy (|7.9p as follows. Choose 
any i G [z^] such that vr'" (z) is undefined and set vr'(j) = i. By (EGC) there is a unique 
i' with 

(-f^s",n")j,* ° {HT.",n")v+i,* = (-f^s",n")«',*- (7-10) 

Also, i' is not in [z/] since by (I7.10p {Hy.» jiii)i,u+i{HY," jiii)u-^i,u+i = {Hy." ,n")i' ■, i-' + 1, and 
the left-hand-side is —1 by Equation (|7.6p . Finally, vr'~ (z') is undefined since no other i 
satisfies (I7.10p . So set 7r'(z^ + j) = i' . 
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Part 3R: Choose a' in £*„ with o"'(l) = 1 and (t'([i/]) = [u] so that, for j, ^ G [u\ 

iHi:,n)^j = (^E,n)^+£j and iH^,u)e,^+j = -(^s,n)^+^,^+j-, (7.11) 

where S denotes the composition of first T," then a'. 

How to choose a': This is symmetric to how we chose vr'. 
Now, since i7'([i^]) = [v], Equation 17.81 imphes 

(^s",n)<,(<>)j = (^E",n)^(^),^+j and (^E",n)<,(^+^)j = -(^s",n)^(^+^),^+j, 

or equivalently 

(^s,n),,, = (^s,n),,,+,. and (^s,n),+,,, = -(^s,n),+,,,+,- (7.12) 

By Equations (j7.11|) and (j7.12|) we can take H' to be the first v rows and columns of 
Ht.,u- d 

Lemma 7.6. Let H be a normalised symmetric n x n Hadamard matrix with n > 2 that 
has an entry —1 on the diagonal and satisfies (GC). Then there is a permutation S in 
Sn with S(l) = 1 and a normalised symmetric Hadamard matrix H' satisfying (GC) such 
that -f^EjE = H2 (B H' . T, and H' can he constructed in polynomial time. 

Proof. In the proof of Lemma 17.51 note that we can ensure 11 = S. If Ha^a = ~1 then 
i = j = 1 in Part 1. □ 



Define H^ as follows. 



Ha 



/+ + + +\ 

+ + - - 
+ - + - 

v+ - - +/ 



Lemma 7.7. Let H he a normalised symmetric n x n Hadamard matrix with n > 2. 
Suppose that H has a positive diagonal and satisfies (GC). Then there is a permutation 
S G (Sn with S(l) = 1 and a normalised symmetric Hadamard matrix H' satisfying (GC) 
such that -ffs.s = H^ (B H' . Ti and H' can be constructed in polynomial time. 

Proof. By Lemma l2. II we know n is a power of 2, say n = 2 ~*"'^. The lemma is trivial for 

k = since H = H4 and S can be taken to be the identity. So suppose A: > 1. Let u = 2^. 

Part 1: Choose S' in Sn with S'(l) = 1 and T.'{iy + 1) = z^ + 1 so that, for j G [2i/]. 

(i?E',s').+i,, = +1 and (i?s',s').+i,2.+, = -1- (7.13) 

How to choose S': We construct a sequence of permutations ctq, . • • , o'2i/ where (Tq is the 
identity and we let S' = a2u- Let H^ denote i?o-j,o-j. For j G {1, . . . 2v}, we define aj as 
follows. If HI_^_^ 2u+j ~ ~^ then aj = o"j_i. Otherwise, there is an 1 < ^ < 2z^ + 1 with 
i ^ i^ + 1 with Hl~^^ = —1. So aj is the composition of first applying aj-i and then 
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w-^ 

"■u+1,1 



transposing 2u+j and i. To see that such an i. exists, note that H^^i^^, has exactly 2i/ ones. 
However, since H is normahsed and has a positive diagonal, Hu-^i^i = Hyj^i^^j^i = +1 so 

Observation: Since -ffs'.S' is Hadamard, (-f^s',S')2i/+i * ^^^ ^'^ positive entries (since its 
dot product with row 1 is 0). Also, half of these are in the first 2u columns (since its dot 
product with row z^ + 1 is 0) . 

Part 2: Choose u' in 5„ with cr'(l) = 1, a'{v + 1) = i' + 1, a'{2v + 1) = 21^ + 1 and 

cr'([2i/]) = [2u] so that, for j G [u], 



{Hy,",T,") 2y+lj - (^S",S")2i.+l,2i/+j - +1 ^'^'^ (-^S",S")2i.+l, 



!/+j 



(i^: 



i;",s'' 



'2u+l,3iy+j 



-1, 



(7.14) 



where S" is the composition of S' then a' . 

How to choose a': We construct a sequence of permutations a[, . . . , a'2y where a'^ is the 
identity and we let a' = (72^,. Let H^ denote H^'-^/^^ij^r. Note that ^2i/+i i/+i = ~^ by 

()7.13p and symmetry of H^. For j G {2, . . . i^}, we define o"'- as follows. If -f^2;7+i i/+j ~ ~^ 
then o"' = o"'„]^. Otherwise, by the observation at the end of Part 1, there is an 1 < £ < u+l 
= — 1. So o"' is the composition of first applying cr'j_i and then transposing 



with H^^l^ 



j and L For j £ {i/ + 1, . . . 2^}, we define aj as follows. If -^21^+1 



2u+j 



-1 then 



0"' = cr'j-i- Otherwise, by the observation at the end of Part 1, there is an 2i/+l < i < 3iy+l 

with ^2iH-i ^ ""■*■■ ^0*^1' i^ ^^^ composition of first applying (t',^ and then transposing 
2i/ + j and i. (The reason that £ > 2z^ + 1 is that the diagonal is positive.) 
Note that S"(l) = 1. Since a'{u + 1) = z^ + 1 and a'{[2u]) = [2u], 



(-f^S",E")i.+ij - {Ht.',T,')„, 



{^+i),'^'(i) 



{H: 



S',S' 



Vi.+l,a'{i)' 



SO Equation (j7.13p gives us 

VjG[2z.],(/7s",S").+ij 



+1 and (-H's",S");.+i. 



2;/+i 



-1. 



(7.15) 



Equations ()7.14|) and ()7.15|) are summarised by the following picture, which takes into 
account the symmetry of Hy," s" • 



Hy. 



S".S" 



/ + ... + 

+ 


+ ... + 
+ 


+ ... + 
+ 


+ ... + \ 


+ ... + 
+ 


+ ... + 
+ 






+ ... + 
+ 




+ ... + 
+ 




+ 
\ + 






/ 
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Part 3: Choose a" in S^ with a"{l) = 1, a"{u+l) = u+1, a"{2u+l) = 2i^+l, a"{[u]) = 
[ly], a"{{u + 1, . . . , 2i/}) = {iy + l,...,2iy} and a"{{2u + 1, . . . , 3u}) = {2i/ + 1, . . . , 3i/} so 
that, for j G [z^], we have the following, where T, denotes the composition of S" then a" . 

(^E,e),,, o (i:fs,E)2,+,. , = (^E,S)2,+1,, (7.16) 

(-f^S,s)^+j^^ o (Fs,e)3j,+j_, = (i^E,s)2j.+i,* (7.17) 

(ffE,s),. , o (Fs,e),+,. , = (^E,s),+i,, (7.18) 

How to choose a": Note that H satisfies (EGC) by Corollary 17.31 hence Hy^'iji" satisfies 
(EGG) (permuting does not change (EGG)). For j G [i^], do the following. Let ii be the 
smallest element in [v] such that the inverse of ii under a" is still undefined. (For j = 1, 
a" is still completely undefined so we will have ii = 1.) Let ^2 be the solution to 

(-f^s",s")«i,* ° iH-E"^T,")u+i,* = (-f^s",s")i2,*- (7.19) 

This equation implies that 

(-f^S",E")n,!^+l(-f^S",S")i'+l,i^+l = {Hs",T.")i2,u+l 

and 

(-ffE",E")n,2i/+l(-f^S",S")!^+l,2!/+l = (-f^E",S")«2,2i/+l. 

Applying Equations ()7.14p and (j7.15p , the left-hand-side of the first of these equations is 
-|-1 and the left-hand-side of the second of these equations is —1, so ^2 £ {i^ + 1, . . . , 2i^}. 
Also, since no other ii satisfies Equation (j7.19p for this value of 12, the inverse of 12 under 
a" is still undefined (so there is no problem with defining it now). Let ^3 be the solution 
to 

(-f^S",E")ji,* ° (-f^S",E")2i/+l,* = (-f^E",E")j3,*- 

This equation implies that 

{H-E",-E")ii,u+l{HT."^T,")2u+l,u+l = {HY.",Y.")i3,u+l 

and 

(-f^E",E")n,2i^+l(-f^E",E")2!^+l,2!/+l = (-f^E",E")i3,2i^+l. 

Applying Equations (I7.14p and ()7.15p . the left-hand-side of the first of these equations is 
— 1 and the left-hand-side of the second of these equations is +1, so i^ G {2v + 1, . . . , 3z^} 
and the inverse of is under a" is still undefined. Let Z4 be the solution to 

(-f^E",E")j2,* ° {Ht." ,T.")2u+l,* = (-f^E",E")i4,*- 

This equation implies that 

(-f^E",E")i2,;^+l(-f^E",E")2i/+l,i/+l = {Hy.",Y,")u,u+1 
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and 

(-f^S",i;")«2,2i'+l(-f^S",S")2i^+l,2jy+l = (-f^S",E")i4,2i/+l- 

Applying Equations (j7.14p and (|7.15p . the left-hand-side of the first of these equations is 
—1 and the left-hand-side of the second of these equations is —1, so 24 G {3i^-|-l, . . . , 4i/}and 
the inverse of 14 under a' is still undefined. Let cr"{j) = ii, a'^u + j) = Z2) (y"{2v + j) = i^ 
and (t"(3i/ + j) = 14^. Note that the choices of ii, 12, 13 and ^4 imply the following, which 
imply Equations (fTTB . (fTTT]) and (TTlHIl . 

(^S",S")a"(j),* ° (^S",S")a"(2i.+i),* = (^S",S")^"(2,.+1),* (7-20) 

(^S",S")<,"(^+j),, ° (^S",S")<,"(3^+i),* = (^S",E")^"(2^+i),* (7.21) 

(^s",s")^"(j),* ° (-H's,s)<^„(^+j)^^ = (^s",E")^"(j.+i),* (7.22) 

Since a'^v -M) = z^ + 1, o-"(2z^ + I) = 2z^ + 1, a"(M) = [i/], a"{{i^ + l,...,2z^}) = 
{1/ + 1, . . . , 2z^} and a"{{2u + I, . . . , 3i/}) = {2z^ -M, . . . , 3u}, Equations (mill and (1713]) 
give us 

Vj G [i^], (^e,e)2^+ij = (^E,s)2^+i,2^+j- = +1 and (H^,^)^^^^^^^^ = (^2,2)2^+1,3^+^ = "1' 

Vj G [21.], (i/s,E),+i,,- = +1 and (Fs,s"),+i,2.+i = '^^ ° 

These, together with Equations (17.160 . ()7.17p . and ()7.18p and the symmetry of -ff^.s, give 
us the result, where H' is the first v rows and columns of -ffE.s- 

Lemma 7.8. Let H be a normalised Hadamard matrix of order n = 2 which satisfies 
(GC). Let X = {Xi, . . . , Xk), Y = {Yi, . . . , 1^). There are index mappings p^ : F^' —>■ [n] 
and p'-' : ¥2 —>■ [n] with p^{0, . . . , 0) = p^'{0, . . . , 0) = 1 and a permutation i: £ Sk such 
that H is represented by the polynomial X^^Y . If H is symmetric then p = p . p , p 
and vr can be constructed in polynomial time. 

Proof. The proof is by induction on k. The base case is A; = 1 for which H = if 2- In this 
case, we take the index mapping p given by /> (0) = 1 and p {!) = 2. p = p and tt is 
the identity. 

For the inductive step, first suppose that H is not symmetric. By Lemma 17.51 there are 
permutations S,n G 5"^ with 5](1) = 1 and n(l) = 1 and a normalised Hadamard matrix 
H' satisfying (GC) such that ifs,n = if 2 © H' . These are constructed in polynomial time. 
By induction, we can construct index mappings p§_i '■ F2~ —>■ [2^~^] and p'^_i : F2~ —>■ 
[2''"^] with /of_i(0, . . . , 0) = /Ofc'_i(0, . . . , 0) = 1 and a permutation vr' G Sk-i such that H' 
is represented by the polynomial 
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Now take />^(Xi, . . . ,Xfc) = ^{I'^-^Xk + p^.iX,, ■ ■ ■ ,Xk-i)) and p^iYi, . . . ,Yk) = 
Il{2''~^Yk + p'j^_i{Yi, . . . , lfc_i)) and let vr G S'fc that the permutation that maps k to itself 
and applies vr' to 1, . . . , A: — 1. 

Next, suppose that H is symmetric and that it has an entry —1 on the diagonal. Using 
Lemma 17.61 we proceed exactly as before except that we are guaranteed (by Lemma l7.6p 
that n = S and that H' is symmetric. Thus, by induction, we are guaranteed that 
p'^_i = Pk-i- So the construction above gives p^ = p^. 

Finally, suppose that H is symmetric and that it has a positive diagonal. Note that 
n > 2. By Lemma 17.71 there is a permutations S G S'„ with S(l) = 1 and a normalised 
symmetric Hadamard matrix H' satisfying (GC) such that i^s,n = H/^ © H' . These 
are constructed in polynomial time. By induction, we can construct an index mapping 
p' : ¥2~ -^ [n] with p'(0, ...,0) = 1 and a permutation vr' G Sk^2 such that H' is 
represented by the polynomial 

Now take p(Xi, ...,Xk) = ^i2''~^Xk + 2'''^Xk-i + p'iXi, . . .,Xk^2)) and let tt e Sk 
that the permutation that transposes k and k — 1 and applies vr' to 1, . . . ,k — 2. 

Proof (Proof of Lemma \2.4\ the Polynomial Representation Lemma). Let n = 2 . Since 
H is positive for A-^ and A*-^, choose a and b such that Ha^ = +1 and (1) a G A-^ or 
A^ = 0, (2) 6 G A*^ or A^ = 0, and (3) If H is symmetric and A^ = A^ then a = b. 
Now let S be the transposition (l,a) and let H be the transposition (1,6). Not that 
(-f^s,n)i,i = +1- Let H be the matrix defined by 

Note that H is normalised. Also, it is Hadamard and it satisfies (GC) since H^jj is 
Hadamard and satisfies (GC). 

By Lemma 17.81 we can construct p^, 'pP and vr such that H is represented by the poly- 
nomial h{X,Y^ := X-^Y . By the definition of "represents", we have 



HpR[-^),pC(y) = -1 ^^ /i(x,y) = 1. 

)p«{x),l 



Define ^^(x) = 1 if (-ffs,n)^fl/-„i i = —1 and ^^(x) = otherwise. Define g^ {y) = 1 if 



(-ffs,n)i pCfy) = —1 and g (y) = otherwise. Now, note that 

(^s,n)p«(x),pc-(y) = -1 ^^ /i(x,y) ©5^(x) ©5^(y) = 1. 

Now let /o^(x) = S(p^(x)) and let /o'^(y) = H(/?-^(y)). Note that H is represented by 
/i(x, y) © 5^(x) © g'~^ {y) with respect to p^ and p^ . 

From LemmaEH p^(0, . . . , 0) = 1 so p^(0, . . . , 0) = a. So if A^ / then /)^(0, . . . , 0) G 
A^. Similarly, /9^(1, . . . , 1) = 6 so if A'^ / then p^{Q, . . . ,^) G A<^. 

Finally, if H is symmetric then Hy,ji is symmetric so H is symmetric so Lemma 17.81 

guarantees that p = 'ff^ . Thus, if A = A , then a = b so T, = Yi so g = g and 
pR ^ pC _ 
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7.4 Linearity 

Proof (Proof of Lemma \2.5[ the Linearity Lemma). Let H he an n x n Hadamard matrix 
and A , A C [n] subsets of indices. Let M, A be the bipartisation of H, A and A and 
let m = 2n. Suppose that (GC) and (R) are satisfied. Let n = 2^ hy Lemma \TJ\ We 
will construct a matrix C and and a reduction EVAL(C, /m,-fm;A) ^ EVAL(M, /m)-fm;A)- 
We will show that EVAL(C, 1^, /m;A) is #P-hard unless (L) is satisfied. 

The reduction is as follows. Let G = {V, E) be an input to EVAL(C, Lm,Lm]A)- We con- 
struct an input G' to EVAL(M, Lm,Lm-A) as follows. Each edge {u, v} G E corresponds to a 
gadget in G' on vertex set {u, v, w, w' , w"} and edge set {{u, w}, {v, w}, {w, w'}, {w' , w"}}, 
where w, w" and w" are new vertices. 

Now let us construct the matrix C. Let F denote the graph with vertices u and v 
and a single edge between them. Clearly, Ga,b is equal to the contribution to Zm{^') 
corresponding to those configurations ^ with ^{u) = a and ^{v) = b. Thus, if c, d, and e 
denote the choice of spins for vertices w, w' and w" , respectively we get 

m mm 

Ga,b = Yl Ma,cMi,^c{Lm;K)c,c Y^ ^ Mc,dMd,e{Im;A)e,e. (7.23) 

c=l d=l e=l 

Here we use that the vertices w, w" have odd degree and the vertex w' has even degree. 
From the definition of bipartisation, we find that Ga^h = Gh^a = for all a E [n] and 
b £ {n + 1, . . . , 2n}. Furthermore, for a,b G [n], 

n n n 

C-a,b = / ^ Ma^n+cMi)^n+c{Im;A)n+c,n+c / ^ / ^ ^n+c,dMd^n+e^m;AJn+e,n+e 
c=l d=l e=l 

n 

= 2_^ HgcHb^c / ^Hd,cHd,e- 
c.eGAC d=\ 

Now, by (R), there are bijective index mappings p^ : Fg — > [n] and p'^ : F2 — > [n] and 
a permutation vr S 5^ such that (w.r.t. p and p ) the matrix H is represented by the 
polynomial h{X, Y) = X^^Y © g^{X) © g^{Y). Let t^ be the inverse of p^ and r*^ be the 
inverse of p'^ . Let L*^ = t'^{K^) and L^ = r^(A^). Also, let a-^ = T^{a), (3^ = T^{b), 
^c = T^ic), 6^ = T^{d) and e^ = r'^(e). Thus, 



Similarly, we get 



HdcHde = (-l)9''(^'')®3''(^'')®^^-(^''®^''). 
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So, for a,b £ [n] 



-y'^ ■{a^®l3^)®9'^ il^)®g'^ {e'=^) V^ /'_i \<5?-{7^®£^) 



Now note that 



, otherwise, 



so for a, 6 G [n], 

(7.24) 
Similarly, 

n n n 

Ca+n,b+n = / ^Ma+n,cMh+n,c{Im;k)c,c / ^ / ^^c,d+nMd+n,e{Im;k)e,e 
c=l d=l e=l 

n 

= 2^ Hc^aHc^b 2_^ Hc^dHe,d, 
c,eGA« d=l 

SO taking a*^ = T'^(a), Z?*^ = T^{h), and 7^ = t^(c), we get 

7«eL« 

Let A^ = \L^\ and A^ = l-L-^]. We will now assume that EVAL(C, /„, ^m;A) is not 
T^P-hard. Using this assumption, we will show that L and L are linear subspaces of Fg, 
which implies that (L) is satisfied. We give the argument for L*-". The argument for L^ 
is symmetric. 

If L is empty then it is a linear subspace of Fg, so assume that it is non-empty. 
Condition (R) guarantees that, since A*^ is non-empty, p*^(0) S A*-^. Hence, G L*-" . 

Let C be the subspace of F2 spanned by L . C contains all linear combinations of 
elements of L . We will show that L = C, so L is a linear subspace of Fg. 

By Equation (j7.23p . the matrix C is symmetric. By Equation (j7.24p . we have Ca,a = nX^ 
for a S [n]. Thus, by Lemma WM (due to Bulatov and Grohe) Ca^b S {— nA'-^,0, nA^} for 
all a, 6 G [n]. Otherwise, EVAL(C,Im,/m;A) is #P-hard. Let x = «? © /3^- Since 
Ca,h G {— nA'-^,0,nA'^}, Equation ()7.24p . implies that for every such x £ "^2^ 

E (-ir^e{-A^,0,A^}. 
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) Hi UilC OUillillclljlUil 1& ^ — ±J ^ ±j &U UilC UUH^UlilC — / ^ 

not possible. Therefore, we get 



Since G L , one of the items in the summation is (— l)'''^ = 1, so the outcome —A is 



Y, i-^V'' e {0, A^}, for all x G ^l (7.26) 

Let Ho = {x £ 1^2 I ^7 G L^\X'7 = 0}- If X ^ '=■0 then x-7 = for all 7 G £. Otherwise, 
by the linearity of C, 

\{^ G C : x-1 = 0}\ = \{-f £ C : x-7 = l}\- 
Thus 






m , if X G So 
, otherwise 

7G£ 



Hence (the characteristic functions of) the sets L'-'' and C have the same Fourier transform, 
up to scaling. It follows that C = L and L is a linear subspace of F2 as required. 

Finally, note that it is easy, in polynomial time, given H, to construct C and to determine 
whether, for all a, 6 G [n], Ca,b £ {— nA'^,0,nA'^} and Cn+a,n+b G {— nA^,0, nA^}. Thus, 
it is easy, in polynomial time, to determine whether (L) holds. □ 

The following fact about linear maps will be useful later. 

Lemma 7.9. Let (j) : ¥2 ^ F2 be a linear map. There is a surjective map / : F2 ^ F2 
and a constant z G N such that 

• /(ci,...,Cfc) • {xi,...,Xi) = (ci,...,Cfc) ■4>{xi,...,Xi), and 

• y{c[,...,c'f,),z = |{(ci,...,Cfc) I /(ci,...,Cfc) = (ci,...,4)}|. 

Proof. Let B be the i x k matrix defining </>, i.e. 0(xi, . . . , xi) = {xi, . . . , xe)B. Define / 
by /(ci, • • • , Cfc) = (ci, . . . , Ck)B'^. Then letting x denote the row vector (xi, . . . , x^), 

f{ci,...,Ck) • {xi,...,Xi) = /(ci,...,Cfc)x^ 

= {ci,...,Ck)B'^x^ 

= (ci,...,Cfc)(x5) 

= (ci,... ,Ck) ■ <f>{xi,... ,xe). 
Fix any c' G F2 and any c G F2 such that /(c) = c'. Note that 

f-\c') = {c + X I X G F^', /(c + x) = c'} . 

As /is linear, we have/ (c+x) = /(c)+/(x) = c'+/(x) so /-^c') = {c + x| xGF^, /(x) = 0}. 
Thus, wetake z = |{xGF^' I /(x) = 0}|. D 
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7.5 The Degree Condition 

Let X = (Xi, . . . , Xk). Every polynomial in g{Xi^ . . . , Xk) £ F2[Xi, . . . , Xk] can be written 
as a sum of distinct monomials of the form Xj^ • Xi^ ■ ■ ■ Xi. for 1 < ii < . . . < ij < k. 
Given a polynomial g{X), let #{g{X)) = \{a G F^ | g{a) = l]\. For a,/3,7 G F^, let 

5a,/3,7(^) = 5(a ®X)®g{l3®X)®i-X. 

Lemma 7.10. Let g G F2[Xi, . . . ,Xj^] be of degree at least 3. Suppose that variables Xy, 
Xs and Xt are contained in a monomial of degree at least 3. Let (3 = and let a G F2 be 
the vector which is all zero except at index r. Then there are polynomials h, /ij.,s, hr^t cmd 
hr such that h is not identically and 

ga,(3,'yiX) = X,Xth{X\{r, s, t})(BXshr4X\{r, s, t})®Xthr,t{X\{r, s, t})(Bhr{X\{r, s, t})©7-X 

(7.27) 
for all J G F2. 

Proof. Let Z denote the tuple X \ {r,s,t}. Let h'{X) be the sum of all monomials of g 
that contain X^, Xg and Xt. Let h{Z) be the polynomial satisfying h'{X) = XrXsXth{Z). 
Note that h{Z) is not identically zero. Choose /ij.,s, /ir,t; hg^t-, hr, kg, ht and /i0 so that 

g{X) = XrXsXth{Z) XrXghr,s{Z) © X,Xt/i,,t(Z) © X,Xt/l,,t(Z) 
®Xrhr{Z) © Xshs{Z) © Xtht{Z) © /lo(^) 

Then for a and (3 as defined in the statement of the lemma, we have 

g{a ® X) ® g{f5 ® X) = g{a®X)®g{X) 

= {{Xr © l)XsXt © XrXsXt) h{Z) 
®{{Xr © l)Xs © XrXs)hr^s{Z) 
®{{Xr © l)Xt © XrXt)hr^t{Z) © hr{Z) 

= XsXth{Z) © Xshr,s{Z) © Xthr,t{Z) © /l^(Z), 



which finishes the proof. □ 

Lemma 7.11. Let g{X) G F2[Xi, . . . ,Xfc]. The following are equivalent. 

1. g has degree at most 2. 

2. For all a and (3 in ¥2, 

• there is exactly one 7 G F2 such that #{ga,i3,'y{X)) G {0,2'^}, and 

. For all i / 7, #(<7a,A7'(^)) = 2'-'- 

Also, if g has degree greater than 2 then there are a and [3 in ¥\ for which there is no 
7 G F^ such that #{ga,f^,^{X)) G {0,2'=}. 
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Proof. Suppose that 5 has degree at most 2. Let g'{X) := g{a(BX)(Bg{(3(BX). Consider any 
degree-2 term XrXs in g. In g' , this term becomes {Xr(Bar){Xs(Bas)(B{Xr(B(3r){Xs(BPs)- 
Now{Xr © ar){Xs ® Os) = XrXg © Xj-Os © OrXg © a^ctsi SO the term X^-Xg cancels in g' . 
We conclude that g'{X) is linear in Xi, . . . , X^ and (2) holds. 

Conversely, suppose that g has degree at least 3. Suppose that variables Xr, Xg and Xf 
are contained in a monomial of degree at least 3. Let /? = and let a € Fj be the vector 
which is all zero except at index r. By Lemma [7. 101 there are polynomials h, hr,s, hr^t and 
hr such that h is not identically and 

5a,/3,7(^) = XsXth{X\{r, s, t})eXshr,siX\{r, s, t})eXthr,t{X\{r, s, t})ehr{X\{r, s, t})©7-X 

Since h is not identically 0, the term XsXfh{X \ {r,s,t}) does not cancel for any choice 
of 7. Hence, there is no 7 such that i^{ga,f3,-f{^)) ^ {0,2^^}, so (2) does not hold. □ 



Lemma 7.12. Let g G F2[Xi, . . . ,Xfc] be of degree at least 3. Suppose that variables Xr, 

1 



Xg and Xt are contained in a monomial of degree at least 3. Let (3 = and let a S Fg be 



the vector which is all zero except at index r. Then there is a 7 S F2 such that 

#{ga,f3,y{X))^2'^-\ 

Proof. Suppose, for contradiction, that i^{ga,i3,-y{X)) = 2^~^ for every 7 G F2. Let Z 
denote the tuple X \ {r,s,t}. By Lemma |7.10| there are polynomials h, hr^s, h^^t and hr 
such that h is not identically and 



ga,PniX) = XsXth{Z) © Xshr,s{Z) © Xthr,t{Z) © /l,(Z) © 7 • X. 



(7.28) 



Let Y £ "^2 denote the vector obtained from 7 by deleting component 7^. Let 7' E 
Fg^ denote the vector obtained from 7 by deleting components 7,., 7^ and 7^. Let 

g'^^p^^iX \ {Xr}) = XsXth{Z) © Xshr,s{Z) © Xthr,tiZ) © /l,(Z) © 7'' " (^ \ {^r}), 

SO ga,i3,'y{X) = g'^^f^^^{X\{Xr})®^rXr. The polynomial g'^f^.^{X\{Xr}) can be simplified 
as in the following table, depending on the possible values of Xg and Xf. 



Xs 


Xt 


9'aS..iX \ {^r}) 




1 



1 





1 
1 


hr{Z)®i -Z 

hr^s{Z)®hr{Z)®i ■ Z®-fs 
hr,t{Z) ® hr{Z) ® j' ■ Z ® -ft 
HZ) © hr,siZ) © hr,tiZ) © hriZ) © 7' 


•^©7s 


©7t 
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Define 



m+ 

Vi- 

V2- 
V3+ 
V3~ 



#{hr{z)ej'-z) 

#{hr,s{Z)(Bhr{Z)(B7' -Z) 

#{hr,s{z) ® hr{z) (B -r' ■ z a 
#ihr,t{z) e hr{z) e j' ■ z) 

#{hr,t{Z) (B hr{Z) (B j' ■ Z (i 



1) 



#(/l(Z) e hr,s{Z) © hr,t{Z) © hr{Z) © 7' • Z) 
#(/l(Z) © hr,s{Z) © hr,t{Z) © hr{Z) © 7' • Z © 1) 



We can express #{g'a,p^^{X \ {Xr})) in terms of %, ?/i+, ry^-, %+, %-, %+ and r/3- 
dep ending on the values of 7^ and 7^. 



Is 



It *{g'^,p,,{x\{Xr])) 





1 

1 





1 
1 



% + r]i+ + V2+ + V3+ 
m + rii' + V2+ + V3- 
% + Vi+ + V2- + V3- 

m + vi- + ^2- + ^3+ 



Since Z is a tuple of A; — 3 variables, each of r/o, r/i+, r/i-, %+, %-, ^^3+ and r/3- is 
between and 2^~^ . Furthermore, we have r/j+ + r/j- = 2^~^ for all i G [3]. We are 
assuming #{ga,i3,'y{X)) = 2^^^ for any 7, so for any 7 with 7^ = 0, Equation (|7.28|) 
implies i^id'a 13 ji^ \ {^r})) = 2^^"^. Altogether, we obtain the following system of linear 
equations, which is applicable for any 7 with 7^ = 0. 



/ 1 1 \ 




/ ^0 \ 




( 2^-^ \ 


110 




' m^ 




2fc-3 


11 




m- 




2fe-3 


110 10 10 




V2+ 


= 


2fc-2 


10 110 1 




V2- 




2fc-2 


110 10 1 




, ^3+ 




2fc-2 


\ 1 1 1 1 / 




\ %- / 




^ 2^=-^ j 



Solving this system yields r/o = ?/i+ = . . . = r/3- = 2^^ ^. Since rjQ = 2^^ ^, 

Vy G F^-3, #(/i^(Z) © 7' . Z) = #(/i,(Z) © 1 © 7' • Z) = 2'=-^ 



(7.29) 



We will use Equation (j7.29p to derive a contradiction. Let W = (Wi, . . . ,Wk-3) and 
let Y = {Yi,...,Yk-s). Let f{W,Y) = hr{W) © hr{Y) ®W -Y. For 7' G F^-^ let 
fy{Y) = f{i,Y). by Equation ^M), 



yj' e ¥',-', #{fy{Y)) 



^k-i 



(7.30) 
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Note that / represents a symmetric Hadamard matrix Hf of order 2 . So equa- 
tion (j7.30p says that all rows of Hf have sum 0. This is impossible because the rows, 
together with the all-ones vector would then be an 2 + 1 element basis of a 2 di- 
mensional vector space. So we have a contradiction. □ 

Corollary 7.13. Let g{X) G F2[Xi, . . . ,Xk]. The following are equivalent. 

1. g has degree at most 2. 

2. For all a ^ 13 in F^, 

there is at most one 7 G F2 such that #{ga,i3,-y{X)) G {0,2 }, and (7-31) 



for all i / 7, #(5„,^,y(X)) = 2'=-!. (7.32) 

Proof. If g has degree at most 2 then (2) holds by Lemma 17.111 Suppose that g has 
degree at least 3. Lemma 17.111 provides an a and /? such that there is no 7 such that 
i^{ga,i3,-y{X)) G {0,2''}. So to prove the theorem we just have to rule out the case that 
every 7 satisfies #{ga,i3,-yiX)) = 2^^^ for this choice of a and /?, and this is ruled out by 
Lemma I7.12[ □ 

Proof (Proof of Lemma \2.6l The Degree Lemma). Let H he, an n x n Hadamard matrix 
and A^, A*^ C [n] subsets of indices. Let M, A be the bipartisation of H, A^ and A*-^ 
and let ra = 2n. Suppose that (GC),(R) and (L) are satisfied. For integers p we will 
construct a matrix C^ and a reduction EVAL(CH) < EYAL{M,LrnJm;A)- We will show 
that if (D) does not hold then there is a p such that EVAL(C'^') is #P-hard. 

The reduction is as follows. Let G = {V, E) be an input to EVAL(C[p]). We construct an 
input G' to EVAL(M, Im,/m;A) as follows. Each edge {u,v} G E corresponds to a "lotus" 
gadget in G' . The vertex set of the gadget is {u,v,u[,v[,u'l,v'l,Xi,yi,z,w \ i G [p]}. See 
Figure [2] for an illustration of the lotus gadget for p = 1. The gadget has the following 
edges, for all i G [p]: {z,Xi}, {w,Xi}, {z,yi}, {w,yi}, {u^u'^}, {u[,u'-}, {xi,ii-}, {xi,vl}, 
{v,v'i}, {Vi,v'/}, {yi,u'i}, and {yi,v'i}. 

Note that the vertices of the gadget have the following degrees: 

d«) = d{v'/) = 1 

d{u'i) = d{vi) = d{xi) = d{yi)) = 4, 

d{z) = d{w) = 2p. 

Furthermore, for the "boundary" vertices u, v we have 

dc'iu) =p-dG{u), dc'iv) =p-dGiv). 

We will stipulate that p is even. Then the degree of vertices, except for the u'^ and v'/ , is 
even. 
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Figure 2: The Lotus gadget for p = 1. 

Now let us construct the matrix C^P' . Let F denote the graph with vertices u and v and 
a single edge between them. Clearly, C^^^ is equal to the contribution to ^M,/m,/m;A(F') 
corresponding to those configurations ^ with ^(u) = a and ^{v) = b. 

By (R), there are bijective index mappings p^ : ¥2 -^ [n] and p*^ : F2 ^ [n] and a 
permutation vr G S'fc such that (w.r.t. p and p ) the matrix H is represented by the 
polynomial h{X, Y) = Xt,Y ® g^{X) © 5''-^(y). Let r^ be the inverse of p^ and r*^ be the 
inverse of p . Let L = r (A ) and L = r (A ). By condition (L) we know that the 
sizes of L^ and L^ are powers of 2. Let |L^| = 2^ and let IL*-^! =2^ .If A^ is nonempty 
then let (f)^ : F2 ^ F2 be a coordinatisation of A^ with respect to p^. Similarly, if A*-" is 
nonempty, let 4> be a coordinatisation of A with respect to p . Let (f) = (j) if A = A 
and this is nonempty and H is symmetric. Note that if A*-" and A^ are empty then (D) 
is satisfied. 

Let Fj be the subgraph of F' induced by {u,Xi,yi,u'^,u'-}. For a,j,6 G F2, let a^ = 
p^{a), c^ = p^i'y) and d^ = p^{6). Let Z^{a, 7, 5) denote the contribution to ^M,/m,-fm;A(Fj 
corresponding to those configurations ^ with ^(n) = a^ and C(xj) = c^ and C{yi) = d^, 
ignoring contributions due to Im;A for vertices u, Xi, and t/j. (We ignore these contributions 
because these vertices will have even degree in G' so these contributions will cancel when 
we use Z{a, (3,^).) Using n + a' to denote the spin at u'^ (which must be in the range 
{n + 1 , . . . , 2n} , otherwise the contribution is zero) and a" to denote the spin at u'- (which 
must be in [n]), we get 

n n 

Z^{a,^,5) = X] X^ M^R,n+a'Mn+a',a"M^R,n+a'MdR,n+a'iIm;A)a",a" 
a'=l a"=l 
n 

a"^ARa'=l 

Plugging in the representation of H where p^{4>^{p)) is the spin a" G A^, we get the 
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following. 



2■^('Q,^^^^^ = ('_]^^S''(")®9''(7)®9'*(<5) y^ (_^_iy'^{'t>'^M) y^ ('_]^AQ'-(a^®0'^{/i),©7vr®<5^)^ 



Note that 



MGFf a'eF| 



a'eF| 



, otherwise 



Equivalently, 






a'GF^ 



n , if (/)'^(//) = a©7©5 
, otherwise 



Thus, Z^{a, 7, (5) =0 unless a © 7 © 5 G L^ and in this case, 

Z^(a,7,5) =n(-l)3''(°)®s«(^)®s«(^)®5''("®^®^). (7.33) 

Our strategy for the rest of the proof is the following: The goal is to prove that either 
there is a p such that EVAL(C[^]) is #P-hard or the following two conditions are satisfied: 

Row Condition Either A is empty or the polynomial g o (p^ has degree at most 2. 

Column Condition Either A is empty or the polynomial g o (p^ has degree at most 2. 

Let us turn to the Row Condition first. Suppose that A^ is nonempty; otherwise there 
is nothing to prove. Let a, 6 G A . Define a and /? in F2 so that (/) {a ) = t (a) 
and c/)^{P^) = T^{b). Note that contribution to ^Af,/„,/„;A(r') of a configuration ^ with 
^{u) = a and ^{v) = 6 is zero unless the spins of vertices u'-, v'-, z and w are in {n+1, . . . , 2n} 
and the rest of the spins are in [n]. Then taking p^ {e) + n as the spin of z and p'^(C) + ^ 
as the spin of tt;, we get 

= E ( E ^''(</'^(«^),7,'5)^^(0^(/?^),7,5)(-l)(^-®^-)-(^®« 

£,CGF§ \7,5gF^ 

Prom Equation (j7.33p we find that if we take any 7' and 5' such that 7' © 5' = 7 © 5 then 
Z^{a,j,S)Z^{P,j,6) = Z^{a,j',S')Z^{P,j',6') for any a and (3. Thus, we can simplify 
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the expression using ip to denote e © (" and rj to denote 7 © (5. 



i/'eFo \7.5eF'= 



2 \/>"^^2 

p 



^ E h E ^''('A''(a''),^,O)Z«(0^(/3^),ry,O)(-l)''- 






^P+1 ^ ^ Z^(<A^(a«),r/,O)Z«(0^(/3^),r?,O)(-l)''- 






12 \7^'i2 



Now, by equation (j7.33p . the contribution for a given 77 is unless (j)^{a^)(Br] and (/)^(/3-^)©ry 
are in L . But (f) {a ) and (/> (/3 ) are in L , so by (L), the contribution for a given 77 
is nonzero exactly when rj S L^. Thus, we can use equation (j7.33p to simplify, writing 77 
as 4>^{fi). 

(J[P] ^ ^P+l y^ [ y^ ^2/_-,_x(;«(9i«(a«))e<;«(</'«(/3«))®9^(</'-^("^)er?)®9«(?i«(/3«)ffir,)er7^-^ ^ 



n3P+l((-l)9 (-^ (" ))®9 ("^ (^ ))) y" y" (-1) 



V^eF^ \^gF|« 



^2 



xg«(<^«(a«))©S«(0«(/3«))V _ , -r. . „ ., ^ l^.,oa.i+,. ^f ^R 



Since p is even, we have I (—1)^ IP ^M = ]^^ Using the linearity of and 

inverting vr, we further simplify as follows. 

^3p+l y^ I y^ /_2\S^</'-"("^ffiM)ffi9^0-^(/9^©/i)©0-"{A')-X 
X6F* \/.GFf 

Since c/)^ is linear, by Lemma 17.91 there is a surjective map / : Fg — > Fg and a constant 
K G N such that </> (/i) • x = /(x) ' ^ &iid for any 7 G Fg the number of x with /(x) = 7 
is K'^ so we can simplify. 

qIp] ^ n^P+'^i^R y^ I y^ /_-^N5^</'-^{a^®M)®9^0-^(/3^©/i)©/i-7 

7eFf UeFf 
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Let 



(7 b] 



(7W 



„3p+l . ^R ■ 

Clearly EVAL(CW) = EVAL(CW). We^will now show that g^ o 0-^ has degree at most 2 
or there is an even p such that FiYAL{C'-P') is ^P-hard. First note that C^^' is symmetric 
and 



c&= E E (-1: 



fi-y 



7GFf \/.eFf 

For X = {Xi, . . . , X^r) and a,h ^ A^ and 7 G Fg , define the polynomial 

For all a,b G A^ we define: 

Ca,b ■■= {7GFf |#a,6,^(X))G{0,2^''}} 

0a,b ■■= {7eFf I #a,5,,W)^ jo, 2^^-1,2^"}} 

W.,, := {7GFf |#a,6,^(X))=2^"-i}, 

where #{ga,b,'y{X)) denotes the number of x G F2 such that 50,6,7(2^) = 1- 

For every 7 G C/^^;, define ^;a,b,7 := Y2 pp£«(~l)^"'*'^ ) which, by definition, satisfy 

Za,b,'y / and \za,b,y\ < 2^ . Let z™^"" = max^gg^,, \za,brf\ and 2™^° = min^eg^,, \za,b,'y\- For 
a,b ^ A , we can simplify the expression for C^^. 



C, 



I \ P 

a,.- E E (-1)'"'"^^'^ 

7eFf \mgf| 



[p] 



Ff 



E2^^'^+ E(-«A7r+ E 

l,7eCa,6 7eGa,6 l&'Ha,b 



|C,,,|2^''^+ E (^^,^7) 
7eGa,f, 

Since p is even, {za^b.-yY i^ positive for all 7 G ^a,6 and thus C^\, is non-negative for all 
a,b ^ A^. If A^ is empty then the relevant condition in (D) is satisfied, so suppose that 



it is nonempty. We will now show that g 



R ^ aR 



even p such that C^^ has a block of rank at least two. 



has degree at most 2 or there exists an 



Case A. There are a, 6 G A-^ such that Qa,b / Choose such a, h. The principal 2x2 
submatrix of C^^l, defined by a and h has determinant 



'-"6,a ^b,b 



r<\p\ ni^p 
'-'a,b ^ 



22^«p _ (^[P])2_ 



(7.34) 



40 



If the determinant is zero, then 



C. 



b] 



2^^p 



1. We consider two cases. If Ca b = 0, then 



C. 



\p\ 



a,b 



\J2-yega,b^^"-^ 



ha) 



2i^P 




^fRp 




ja. 


?>i(<rF 




2£flp 




<2^" 


^ max X p 




U^V 




<2^" 


/2^^-iy 




I 2^^ J 




<2^" 


•e-^/^^". 



£«N 



(because z^ < 2^^ ) 



This is less than one for ah p > i^2^ . Hence the determinant (j7.34p is nonzero. Fur- 
thermore, as Qa,b 7^ we have C^j^ ^ and hence C^^' contains a block of rank at least 
two. This implies the #P-hardness of EVAL(C''p') by Lemma 16.41 (Recall that C^f^ is 
non- negative since a,b (z A .) 

For the other case, suppose \Cab\ ^ 1- Then 



C. 



[p] 

a, 6 



2^-Rp 



-t^p 



> \CaA > 1- 

Here, the second-but-last inequality holds, because z™™ > and (by the precondition of 



C 



[pl 



case A) Qa^ / 0. Hence again we have —^ / 1, and the determinant (I7.34p is nonzero. 
As in the first case this implies the #P-hardness of EVAL(C'l^]). 



Case B. For all a,b £ A^ it holds that Qai 



Then for all a,b £ A^ we have 



C. 



[p] 

a,b 



\C .12^ P+ V z^ -\C u\2^ 
■y&Qa,b 



So the principal 2x2 submatrix of C^^^ defined by a, b has determinant 

_ r,2e.^pf-\ \n \2\ 



/=^[p] ri[p] 

"-"fe,a ^bh 



2i^P 2^''P\Ca,l 
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This determinant is zero if and only if \Ca^b\ = 1; and the submatrix is part of a block iff 
C-a,b 7^ 0- Hence, we have #P-hardness by Lemma 16.41 if there are a,b (z A^ such that 
\Ca,b\ ^ {0, 1}. Assume that for all a, 6 e A^ we have \Ca^b\ G {0, !}• Define sets 

T:={(a,6) | a G A^, 6 G A^, |Ca,6| = l,a/6}, 
Z := { (a, 6) I a E A^, h G A^, JCbl = 0, a ^ b] . 

Obviously, these form a partition of pairs of distinct elements in A . In other words, for all 
a / 6 e A-^ there is at most one 7 G F2 such that ^{gafi^-yiX)) G <^ 0, 2^ >. Furthermore, 

Qa,b = implies that for all other 7' / 7 we have #{ga,b,'y'{X)) = 2^ ^^. But Corollary 
17.131 implies that in this case g o (p^ has degree at most two. This finishes Case B and 
hence the proof of the Row Condition. 

For the Column Condition, in a symmetric way to how we defined Z (a, 7, (5), we let 
Z^{a, 7, S) denote the contribution to •^M,/„,/„;A(ri) corresponding to those configurations 
^ with ^(n) = n + a*-^ , ^{xi) = n + c'~^ and ^{yi) = n + dP , ignoring contributions due to 
Im-.k for vertices u, Xj, and yj. Using this, we can compute C^4_„„^^ for a,b (^ A^ and 
show that, if A*^ is nonempty, then either g'^ (\P has degree at most 2 or EVAL(C[^]) is 
#P-hard. 

Finally, we note that it is straightforward, in polynomial time, to determine whether 
EVAL(CH) is #P-hard or (D) holds. D 

Corollary 7.14. Let H be a symmetric n x n Hadamard matrix and A = A C [n] 
identical subsets of indices. If H is positive for A^ and A^' then EVAL(//, /„, /„.^ii) is 
polynomial time computable if, and only if, H A^ and A'-" satisfy the group condition 
(GC) and conditions (R), (L), and (D). Otherwise EVAL{H, In, In-An) is #P -hard. If 
H is not positive for A^ and A^ then FiYAL{H, I^, In/^^n) is polynomial time computable 
if, and only if, —H A and A satisfy the group condition (GC) and conditions (R), 
(L), and (D). Otherwise FiVAL(II, In, I^-ar) is ^F -hard. 

Proof. By the equivalence of EVAL(//, /„, /^.^u) and EVAL(— i7, /„, In-A_R) we can assume 
that H is positive for A and A . First, suppose that one of the conditions is not satisfied. 
By Theorem[221 EVAL(M, Im,Im;k) is #P-hard. Since M is bipartite, EVAL(M, Im,Im;A) 
remains ^^P-hard when restricted to connected bipartite instances G. But for these in- 
stances, Zm,i^j^.a{G) = 2ZHjr.,i^.!^R{G), so ¥MAh{H,In,In-KR) is #P-hard. 

It remains to give the proof for the tractability part. For symmetric H and A = A 
satisfying (GC), (R), (L) and (D), we shall show how to compute Zh,i„,i r{G) for an 
input graph G in polynomial time. Let Vo ^V denote the set of odd-degree vertices of G 
and Ve = V\Vo. We have 

Zh,i„,i^,^r{G) = 5Z n ^CH.CM n (^";A«)CM,CM = Yl n -^^W.CM 

C.V^[n] {u,v}eE veVo ^-.V^ln] {u,v}eE 

Fix a configuration ^ : V ^>- [n] and let p = p^ = p^ he the index mapping and h the 
F2-polynomial representing H as given in condition (R). Let furthermore (f) := (p = (p 



42 



be the coordinatisation of A as given in condition (D). Let r be the inverse of p and 
L = r(A^). Then ^ induces a configuration (^ : 1/ — > F^ defined by <;" = r o ^ which imphes, 
for ah u,v £ V that /i(<j(u),<j(f)) = 1 iff H^^u),({v) = ~1- We can simphfy 

Zhi I (G) = y^ TT (-l)^(^°?(")'^°'^M) = y^ (_i)®{..,«}6B^('?(").?M) 

C:y^[n] {u,v}eE <j:y^F§ 

€(Vo)CA« ?(K,)CL 

(7.35) 
Define, for each v £V a tuple X'" = {Xf, . . . , X^) and an F2-polynoniial 

{u,v}eE 

Let var(/iG) denote the set of variables in he and, for mappings x '■ var(/iG) — > F2, 
we use the expression x(^") •= (x(-'^r)) • • • > x(-^a:)) as a shorthand. Define /ig(x) := 
©{■ui;|ge ^(x(-'^"))X(-'^'')) and note that this is a sum in F2. 

For a G F2 let 

Sa := \{x ■■ var(/iG) ^ F2 | xiX"") G L for ah 7; G K and hcix) = a}\- (7.36) 

Hence, by equation (j7.35p . Zh,i„,i j, (G) = sq — si. It remains therfore to show how to 
compute the values Sa- Clearly, 

/iG = {X^)nX''(Bg{X^)(Bg{X^)= {X^).X^ (B ^ g{X^) 
{u,v}ge {u,v}eE veVo 

as the term g{X'") occurs exactly deg(f) many times in the above expression and thus 
these terms cancel for all even degree vertices. 

By equation (I7.36P we are interested only in those assignments x which satisfy xi^^) ^ L 
for all I' G Vo- With |A^| = 2' for some appropriate /, we introduce variable vectors 
V" = {Y^, . . . , yj'') for all u G VJ,. li u £ Vo or v G Vo then we can express the term 
{X^)tjX^ in ho in terms of these new variables. In particular, let 

h'6 = {HYnh-HYn® (x«)^.x-© (0(y«))^.x^ 



Let 



{u,v}eE {u,v}gE {u,v}£E 

U,VGVo U,weVe u£Vo,v£Ve 



h'a = h'l,®^®g{ct>{Y^)). 

v&V° 



Then we see that 



Sa := |{x : var(/i'G) ^ F2 | h'cix) = a}\. (7.37) 

By condition (D) g o (j) \s a polynomial of degree at most 2 and therefore K'q is a 
polynomial of degree at most 2. Furthermore, we have expressed Sa as the number of 
solutions to a polynomial equation over F2. Therefore, as in the proof of Theorem 12. 2^ the 
proof now follows by Fact 12.71 □ 
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Proof (Proof of Theorem \l.^) . Let H he a symmetric n x n Hadamard matrix and A = 
A*^ = [n]. Then H is positive for A-^ and A*^. Let M, A be the bipartisation of H, A^, A*^. 
Suppose first that H has no quadratic representation. Then there are no index mapping 
p = p^ = p^ and coordinatisation (^ = (j)^ = (j)^ such that conditions (R) and (D) are 
satisfied. Hence by Theorem 12. 2 1 EVAL(M, /m;-fm;A) is #P-hard. Since M is bipartite, 
EVAL(M, Im,/m;A) remains #P-hard when restricted to connected bipartite instances G. 
But for these instances, ZM,im,im{G) = 2ZH,i„,in{G), so EVAL(ii', /„,/„) is #P-hard. 
Suppose next that H has a quadratic representation with index mapping /) : F2 — > [n] 
and polynomial h{X,Y). Instead of going through Theorem 12.21 it is easier to prove the 
tractability of EVAL(i?) directly along the lines of the proof of the tractability part of the 
theorem. We leave the details to the reader. This is similar to the tractability part of the 
proof of Corollary 17.141 □ 

8 The Proofs for Section [3] 

8.1 Technical Preliminaries 

Lemma 8.1. Let C E M™^™ he a symmetric matrix and let A+ and A^ he diagonal 
mx m matrices. Let D he the component-wise sum D = A+ + A^ and let O = A^ — A^. 
Let A he the tensor product 

Let A he the 2m x 2m matrix such that, for all i G [m] and j G [m], Ajj = Af-, 
Ai^m+j = Am+ij = 0, and Am+i,m+j = A^.. Then 

Zc,D,o{G) = Zaa{G) for all graphs G. 

Proof. It is useful to think of A and A in terms of four m x m tiles as follows. 

. ( G -C\ ^ ^ / A+ 
^=V-C g) a^dA=^ Q ^_ 

We will simplify the expression for Za,^{G) now. Let ^ : ^ — > [2m] be a map such that, 
for some w ^V , ^{w) G [m\. Let V' be the mapping such that for all v &V 

^^ ' ^^ ^ \ , otherwise. 
Then 

11 ■^■)p{u),il)(v) — 11 ■^np(w),'il){w) 11 ■^ip{w),il)(v) 11 ^V){n),i/)(i)) 
{u,v}&E {w,w}^E {w,v}&E {u,v}<^E 

= n ^e(«').c(«') n ~^5("'):?(^) n ^«w.cw □ 

{w,w}^E {w,v}^E {u,v}&E 
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which imphes that 

{u,v}£E {u,v}£E 

where deg{w) denotes the degree of w in G (self-loops add two to this degree). Since 
Uvev ^av),av) = ^a^),aw) Uw^vgV ^e(f),?(f)' ^^ have 

S,:V^[2m] {u,v}eE vl^V 

^■.V^[2m] {u,v}&E wjtveV 

g{ui)G[r?i] 

As this argument can be applied independently to all w; G y we obtain 

■^A,a(G) =2-^11 ^C{«),?(t>) 11 (^C{'^):C(«') + i~^)^^^ ^^■inHM,rriH(w)) 

S,:V~^[m] {u,v}eE weV 

= J2 H ^a^uH n ^a^uM n '^a^uM 

^■.V^[m] {u,v}gE w€V w£V 

deg(ui) even deg{ui) odd 

= Zc,D,oiG). 

Corollary 8.2. Let C be a symmetric m x m m,atrix which contains exclusively blocks of 
rank 1. Let D and O be diagonal mx m matrices. Then the problem EVAL(C, D,0) is 
polynomial time computable. 

Proof. By Lemma 18.11 the problem EVAL(C, Z),0) is polynomial time equivalent to a 
problem EVAL(A, A) with A a matrix consisting of blocks of row rank at most 1. Thus 
the statement of the corollary follows from Lemma 16. 3i □ 

8.1.1 Extended Twin Reduction 

Unfortunately the Twin Reduction Lemma 16.21 does fully satisfy our needs. As we are 
dealing with possible negative rows we will be in a situation, where it is useful to reduce 
matrices even further, namely by collapsing two rows j4j^* and Aj^^ into one if Ai^^^ = itAj^*. 
To achieve this we say that two rows j4j^* and Aj^^, are plus-Tninus-twins (pm-twins for 
short) iff j4j^* = ±^j,*. This induces an equivalence relation on the rows (and by symmetry 
on the columns) of A. Let Ii, . . . Ik be a partition of the row indices of A according to 
this relation. For technical reasons it will be convenient to partition the sets Ii into the 
positive and the negative part. That is for every i E [k] we define a partition (Pi,Ni) of 
Ii such that Pi ^ 9 and for all v, v' e Pi and /u, /i' G Ni we have Ay^^ = A^i^^, A^^^, = A^i^^, 
and ^/y^jf; ^ ^u,** 
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The pm-twin-v esolvent of A is the matrix defined, for all i,j £ [k], by 
T {A)ij := Afj,^i, for some fi £ Pi^v £ Pj. 

This definition is technical and seems to be counter-intuitive, as we are not taking the 
Ni into account. However its motivation will become clear with the following Lemma and 
it is still well-defined, even though possibly Ni =% for some i G [k]. 

As before, we define a mapping r : [m] — > [k] defined by ^u G /t-(^) that is r maps 
/i G [m] to the class Ij it is contained in. Therefore, we have T^{A)^(j\.^ij\ = ^Aij for all 
i,j £ [m]. We call r the pm-twin-resolution mapping of A. Define A^ = A'^i U . . . U N^ and 
P = Pi U . . . U Pfc. Then in particular 

r^{A),ii),rU) = Aj for all {i,j) G (P X P) U (iV X iV) 

'^^(^)r(i),r(i) = -Aj for all {i,j) e{PxN)U{NxP) 

Lemma 8.3 (Extended Twin Reduction Lemma). Let A be a symmetric mxm m,a- 
trix and A a diagonal mxm matrix of vertex weights. Let (Pi,Ni), . . . , {Pk,Nk) be a 
partition of the row indices of A according to the pm-twin-relation. 
Then 

Zaa^G) = Zr±{A),D,oiG) for all graphs G 

where D and O are diagonal k x k matrices defined by 

Di;i = ^ l^y^v + ^ A^,^ and Oi^i = ^ A^^^, - ^ A^,^ for all i £ [k]. 

Proof. Define Jj = Pj and Jk+i = -Ai for all i G [k]. W.l.o.g. we may assume that if 
there is a minimal I G [k] such that Jk+i = then for all j > I we have Jk+j = (this 
can be achieved by appropriate relabelling of the Pj and Ni). Let / := A; + 1 if all J^+j 
are non-empty. Then Ji, . . . , Jk+i-i are the equivalence classes of A according to the 
twin-relation. Therefore, the Twin Reduction Lemma 16.21 implies that for the diagonal 
{k + £ — 1) X (k + i — 1) diagonal matrix A" defined by A" • = J2ueJ '^^^.'^ ^^ have 

■^a,a(G') = ^r(A),A"(G') for all graphs G. 

Let n' := k + I — 1 and note that by the definition of the sets Ji, T(A) is the upper left 
n' X n' submatrix of the 2k x 2k matrix 

that is T{A) = Mr„/ir„/i. Define a, 2k x 2k diagonal matrix A' such that A^ ^ = A"j for all 
i G [n'] and A^ ^ = for all n' < i < 2k. Then 

Zm,A'{G) = ^r(A),A"(G') for ah graphs G. 
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Moreover, by the definition of A", tlie matrix A' satisfies, for all i € [k], 

Now, by Lemma [STTl Zm,A'{G) = Z'j-±(^y^^^ij/Q, where D' and O' are k x k matrices such 
that Dl, = ^[. + A',^, ;,^. and O^ = A^ .'-'a;^^.^^^.. But by Equation dH]), we see that 
D' = D and O' = 0. D 

Lemma 8.4 (Row-Column Negation Lemma). Let C he a symmetric m x m m^atrix 
and D, O diagonal m x m matrices of vertex weights. 

Let i E [m] and define C as the matrix obtained from C by multiplying row and column 
i with —1. Let O' be the matrix obtained from O by negating the diagonal entry O-i^i. Then 

Zc,D,o{G) = Zc',D,0'{G) for all graphs G. 

Proof. Let G = {V, E) be a graph and Vq, Ve the sets of odd (even) degree vertices in V . 
Recall that 

zc,D,o{G) = X] n qh:?(^) n -^cM.cM n ^€{«).«{«) 

^■.V^[m\ {u,v}&E v&Ve v&Vo 

Fix some mapping ^ -.V ^>- [m]. We will prove the Lemma by showing that 

n Q(«).5M n ^C(^):C(^) n ^iii^)^!^^) = n '^5h,c(^') n ■^5(^),5{^) n ^5m,5w 

{u,v}&E V&Ve. V&Vo {u,v}&E DgVe V&Vo 

Define W := i^^{i) and let We := Vg Pi VF and Wo ■.= Vor\W denote the even and odd 
degree vertices in W . By the definition of O' we have 

V&Vo V&Vo 

Furthermore, for all edges {u, v} £ E we have that G^(^u),^{v) = Glrs g/ n if and only if 
either both u,v & W or u,v ^ W. If exactly one of the vertices is in W then C^(^,) ,c(^,) = 
—G'c,-. f. ^. Therefore, if we denote by e{W, V \ W) the number of edges e = {u, v} in G 
such that exactly one vertex is in W, we have 

n ^ku),iiv) = (-i)'(^'^\^) Yi Q(«),cw- 

{u,v}eE {u,v}eE 

To finish the proof it thus suffices to prove that 

e{W,V\W) = \Wo\ (mod 2). 

The proof will be given by induction on the number \W\ of vertices in W. The case that 
W = 9 is trivial. Assume therefore that there is a vertex w £ W and let U := W \ {w}. 
By the induction hypothesis, we have e{U,V \U) = |C/o|(mod 2). If it; has even degree 
then \Wo\ = \Uo\ and w either has an odd number of neighbours both in U and V \U 
or it has an even number of neighbours in both sets. If otherwise w has \Wo\ = 1 + \Uo\ 
and the parity of the number of neighbours oi w in U is opposite to that of the number of 
neighbours in V \U. This finishes the proof. □ 
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8.1.2 Pinning vertices 

In the proof of Lemma [3.1l it will be convenient to "pin" certain vertices of the input graph 
G to prescribed spins. We will develop the tools which are necessary for this now. These 
results extend analogous techniques used in [U] and [B]. 

Let A he an m X m matrix and D a diagonal m x m matrix of positive vertex weights. 
In the following, a labelled graph is a triple G = {V,E,z), where {V,E) is a graph and 
z € V. For a labelled graph G = (V, E, z) and a. k £ [m], we let 

ZA,D{k,G) = {Dk^kT^ 2Z n ^C(«).€W ■ n -^SM'SM 

i(z)=k 

The product GH of two labelled graphs G and H is formed by taking the disjoint union 
of the graphs and then identifying the labelled vertices. Let H'^ denote the product of H 
with itself taken s times. Note that ZA^£){k,GH) = ZA^£){k,G)ZA,Dik,H) for all labelled 
graphs G and H. 

Recall that a twin-free matrix ^ is a matrix such that Ai 7^ Aj for all row indices 
i 7^ j. Furthermore an automorphism of {A,D) is a bijection a : [m] — > [m] such that 
Aij = Aq(j),^(j) and Di^i = Da(i)^a{i) for all i £ [m]. The following lemma follows by a 
result of Lovasz (Lemma 2.4 in |21j). 

Lemma 8.5. Let A £ i^™-x™ ^g twin free, D G ]^"^x"i a diagonal matrix of positive vertex 
weights and i,j G [m]. If for all labelled graphs G we have 

ZA,D{i,G) = Za,d{J,G) 

then there is an automorphism a of{A,D) such that j = a{i). 

We furthermore need some standard result about interpolation, which we use in the 
form as stated in [9] Lemma 3.2: 

Lemma 8.6. Let wi, . . . ,Wr be known distinct non-zero constants. Suppose that we know 
the values fi,---fr such that 

r 

/j = ^ Cjw'j for all i£ [r]. 
i=i 

Then the coefficients ci , . . . , c^ are uniquely determined and can be computed in polynomial 
time. 

Lemma 8.7 (Pinning Lemma). Let A G M^^™ be a symmetric matrix and A G M™^"^ 
a diagonal matrix of positive real entries. Then for every labelled graph G and every 
k G [m], we can compute ZA,A{k,G) in polynomial time using an EVAL(^, A) oracle. 

Proof. Let the matrices B and D be the result of twin-reduction (Lemma l6.2|) when applied 
to A and A. In particular, B is twin-free and ZA,AiG) = Zb,d{G) for all graphs G. 
Therefore, using the oracle, we can compute Zb^d{G) in polynomial time (for input G). 
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Consider a graph G = {y,E) with a labelled vertex z and a particular spin k £ [m]. we 
will show how to compute ZB^oik, G) using an oracle for Zb,d- 

Call spins i,j G [m] equivalent if there is an automorphism a of {B,D) such that 
j = a(i). Partition [m] into equivalence classes Ii,...,Ic according to this definition. 
For every spin j in equivalence class /,, let Cj denote the size of the equivalence class — 
Cj = \Ii\. For every equivalence class i € [c] let ki denote a particular spin ki £ li. 

For any two equivalent spins a and a' we have Zb,d{o., F) = Zb,d{o,' , F) for every graph 
F. Therefore, 

c 

Zb,d{G) = Y,Ck,ZB,D{h,G) (8.2) 

We will now prove the following claim. The result follows by taking S = Uiefcli^*}- 

Claim 1. Given a set 5 of inequivalent spins and a spin k £ S we can compute ZB^oi^, G) 
in polynomial time using an oracle for computing Ylkes'^kZB,D{k,G). 

Proof. The proof is by induction on l^l. The base case IS"! = 1 is straightforward, so 
assume \S\ > 1. We will show how to compute ZB,D{k,G) (for any spin k £ S) using an 
oracle for J2kes'^kZB,D{k,G). Fix distinct spins i and j in S. By Lemma [831 there is a 
labelled graph Gij such that 

ZB,D{i, Gij) ^ Zb,d{3, Gij). (8.3) 

Note that the construction of Gij takes 0(1) time since Gij does not depend on any 
input graph G. Partition S into classes Ji, . . . , Jj such that v, v' £ J^ iff Zb,d(t^-, Gij) = 
Zb,d{v' -.Gij). We will show below how to compute Ylik<^j <^kZB,D{k,G) (for any fi £ [t]) 
using an oracle for Ylik&s^kZB,D{k.,G). Once we've done that, we can finish as follows. 
For a fixed k £ S, suppose k £ Jf^. Note that \J^\ < S since one of spins i and j is not 
in J^. By induction, we can compute ZB,D{k^G) using the newly-constructed oracle to 
compute ZlfcGJp CkZB,D{k, G). 

To finish, we now show how to compute J2ke.J '^kZB,D{k, G) using an oracle for J2keS ^kZB,Dik, G). 
For every n £ [t], let s^ be a spin in J^. Let w^ = Zb,d{s^, Gij). Let 

fr = Y.CkZB,D{k,GGlj) 

kes 
= E ^ CkZBMk,GGlj) 

= ^ ^ CkZB,D{k,G){ZB,D{k,Gij)Y . 

= X] ^/ X] CkZB,D{k,G). 

ne[t] keJfi 

Note that we can compute fr in polynomial time using the oracle. Now by Lemma [8.61 we 
can recover XlfceJ '^kZB,D{k, G) for every fi apart from the one with w^ = (if there is a 
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/i with Wf^ = 0). But we can recover this one, if it exists, by subtraction since 

y^ CkZB,Dik,G) = '^CkZB,D{k,G) -'^'^ CkZB,Dik,G). 

keJfj, fees vj^fikeJi, q 

The following Corollary will be helpful in the proof of Lemma 18.121 

Corollary 8.8. Let C G M™^'" he a symmetric matrix and D,0 e M^^™' diagonal ma- 
trices such that the diagonal of D is positive and that of O non-negative such that D — O 
is non-negative. Then, for every labelled graph G and every k G [m], we can compute 
Zc,D,o{k,G) in polynomial time using an EVAL(C, -D,0) oracle. 

Proof. Let A^ and A^ be diagonal m x m matrices with A^^ = (-Dj,j + Oj,j)/2 and 
Ki = i^v - Oi,i)/2. Let Let 

A=( \ ~^) C. 



Let A be the 2m x 2m matrix such that, for all i G [m] and j G [m], Ajj = A^-, 

Zc,D,o{G) = Za,a{G) for all graphs G. 



Aj,m+j = Am+jj = 0, and Am+i,m+j = Ai,i- Then by LemmalO 



Let I = {i £ [2m] \ Aj^j ^ 0}. Since D -\- O and D — O are non-negative, we have that 
the matrix Ajj has a positive diagonal. By inspection we have 

ZAAiG) = ZAjiAiiiG) for ah graphs G. 

By the Pinning Lemma [8. 71 we can compute the value ZAji,Ai[{k, G) by an algorithm with 
oracle access to EVAL(^//, A//). Now, ZAjj,Aiiik,G) = Zc,D,o{k,G) for every k G [m]. 
This finishes the proof. □ 

8.1.3 Tensor Product Decomposition 

The following technical Lemma which will be used in the proof of Lemma 13. 2i 

Lemma 8.9. Given symmetric r xr matrices A and D and mxm matrices A', D' . Then 

Za(S)A',d®d'{G) = Za,d{G) ■ Za',d'{G) for every graph G. 

Proof. We consider the indices oi A<Si A' and D ® D' as, pairs (i, j) G [r] x \m] such that, 
e.g. 

Let vr : [r] X [m] — > [r] and p : [r] x [m] -^ [m] be the canonical projections i.e. for every 
(i, j) G [r] X [m\ we have '7T{i,j) = i and p{i,j) = j. 
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Thus 

^■.V^lr]xlm\uveE veV 

S,:V~*lr]x[m\uveE veV 

g:y^[r] uv€E vGV 

= Za,d{G) ■ Za',D'{G) 

It is not hard to see that this kind of decomposition can be performed for parity- 
distinguishing partition functions as well, as the following lemma shows. 

Lemma 8.10. Suppose that A' is a symmetric m! x m! m^atrix and D' and O' are diagonal 
m' X m' matrices. Suppose that A" is a symmetric m" x m" matrix and D" and O" are 
diagonal m" x m" matrices. Then, for every graph G, 

ZA'®A",D'®D",0't^O"{G) = Za\D',0'{G) ■ Za",D",0"{G)- 

Proof. Let A = A' ® A", D = D' ® D" and O = O' i^ O". We consider the indices of A, D 
and O as pairs {i,j) G [m'] x [m"] such that, for example, 

Let it' : [m'] x [m"] -^ [m'] and it" : [m'] x [m"] -^ [m"] be the canonical projections i.e. 
for every {i,j) G [m'] x [m"] we have 7r'(i, j) = i and ^"{ijj) = j. 

With V^ C y the set of even degree vertices and Ve = V \Vo we have 

za,d,o{g) = XI n ^5{«),?(^) n -^ew.cw n ^5m.cm 

^■.V^[m']x[m"] {u,v}£E vGVe v£Vo 



With 



and 



veVe veVe veVe 

veVo veVo v&Vo 



n ^?(«).cw - n K'iii(u))y(av)) n ^'7^"(5{«)),^"(c(^)) 

{u,v}eE {u,v}<=E {u,v}eE 
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we therefore have 

Za,D,o{G) = ( 2^ 11 \'{u),tP'{v) 11 ^tP'{v),tP'{v) 11 0^'{v),tP'{v) 

^ij':V-*[rn']{u,v}eE vGVe vGVo 

/ , 11 ^i:"{u),iP"{v) 11 ^i)"{v),iP"(v) 11 Oi/;"{«),i/)"{«) 

^•0":y^[m"] {«,-u}G-B »;gV'e vGVo 

= ^A',D',0'i'^) ■ ^A",D",0"('^) 

Corollary 8.11. Let B' he a symmetric m' x m' block and let D and O he diagonal 
m' X m' matrices. Let B" he a symmetric m" x m" hlock and let D and O be diagonal 
m" X m" matrices. Let D^ = D^' D^" and O^ = O^' (g) O^" and B = B' (g) B" . If 



FiYAL{B",D ,0 ) is polynomial time computable then 

EyAL{B,D^,0^) = EVAL{B',D^',0^'). 
Proof. For every graph G, Lemma |8. lUI gives 

^B'(g)B",D«'(g)D«",0«'(g>0«"(^) = -^B',D«',0«'(^) ' ^B",D«",0«" (^)- 

If EVAL(i?", D , O ) is polynomial time computable then this gives 

EYAL{B,D^,0^) = EYAL{B',D^',0^'). 
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Lemma 8.12. Let B' he an m' x n' block, D and O he diagonal m' x m' matrices and 
D^ and O^ be diagonal n' x n' matrices. Let B" he an m" x n" block, D^ and O^ he 
diagonal m" x m" matrices and D and O be diagonal n" x n" matrices. Let 

n' ( D""' M ,n" ( ^''" M ^n ( D^' ® D^" 

^ = V z^C J -'^^ = ( DC" ) -^^= ^ DC'g^DC" 

and let O and O', O" be constructed from O , O and O , O in the analogous way. Let 
A, A', A" he the connected bipartite matrices with underlying blocks B := B' g) B" , B' and 
B" respectively. 

If EYAL{A" , D" , O") is polynomial time computable and D + O and D — O have only 
non-negative entries then 

EVAL(^, D, O) = EVAL(A', D' , O'). 

Proof. Note that Za.d.o{G) = unless G is bipartite. Therefore we will assume in the 
following that all graphs G are bipartite and that {U, W) is a partition of the vertex set V 
into two independent sets. Assume first that G is connected - the case of non-connected 
graphs will be handled later. Note that j4 is a square matrix of order m + n for ra = m'm" 
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and n = n'n". For diagonal rxr matrices D,0 a set X CV and a configuration S, '■ X ^ [r] 
define 



deg(a;) even deg(a::) odd 



By the above definitions we have, 



Za,d,o{G) = X] n '^i{u),i;iw)^D,o{U^O^D,o(W,tp) 



(■.U^[m+n] {u,w}eE 
4>-W^[m+n] 



And therefore, since G is connected 



Za,d,o{G) = X] n ^({u),i,{w)'^DR,oRiU,O^DC,oc(W,ij) 



^■.U^[m] {u,w}€E 



+ X] n ^i^iw),an)'^DC,oc{U,0'^DR,OR(W^^) 



C-U^[n] {u,w}eE 
ip:W ^[m] 



Define 



ZXD,oiG) ■= Y^ n ^a^),'4>M^DR,oi'iU,0^DC,oc{W,ij) (8.4) 

S,:U^[m\ {u,w}eE 

and 

ZXD,oiG) := Y^ Yl ^i'M,({u)^DC,oc{U,O^DR,oRiW,i^) (8.5) 

S,:U^[n] {u,w}eE 
xp:W ^[m] 

That is 

Za,d,o{G) = ZXd,o{G) + ZXd,o{G) (8.6) 

For matrices A' ,D' , O' and A" , D" , O" we define the analogous expressions {Z*^, ^, q/{G), 
etc.). 

We consider the indices of B' (gi B" as pairs. That is row indices are (i', i") £ [m'] x [m"] 
and column indices become ij',j") G [?^'] x [n"]. 

[B ® B )(i',i")(j',j") = B^iy ■ B^iiy, 

Let p' : [m'] x [m"] ^ [m'], p" : [m'] x [m"] -^ [m"] and 7' : [n'] x [n"] ^ [n'],^" : 
[n'] X [n"] -^ \n"] be the canonical projections. That is for {i',i") £ [m'] x [m"] we have 
p'{i',i") = i', p"{i',i") = i" and for {j',j") G [n'] x [n"] we have j'ifj") = f and 
yij',j") = j" ■ Therefore, for sl\ S^ : U ^ [m] and i{j : W ^ [n] we have 

11 -^SH-V-Cw) = 11 -^p'og{M),7'o^{w) ■ 11 ^p"o^{u),'y"otp{w) 

{u,w}eE {u,w}eE {u,w}eE 
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and 






Hence, we can rewrite equation ([8 



/ 



■'A,D,0 



(G) 



Yl n ^C'(«),t/''M^DB',0«'(^'^')^DC",0C'(^>V'') 



. i':U^[m'] {u,w}eE 
\ijj':W^ln'] 

( 



) 



\ 



Yl n ^C"(«),i/'"{u;)^15«",0«"(^'^'VdC",0C"(^>V'") 

. C":C/->[m"] {u,«i}G£; 
\^":VK-^[n"] y 

With an analogous argument this extends to Z'^y) o(^)- ^^ therefore have 

^A£',o(G') = ■^A',D',0'(^) ■ ^A'',D",0"(^) (^•''') 

Claim 1. The values Z^£,q(G) and Z^£)q(G) can be computed in polynomial time for 
every graph G by an algorithm with oracle access to EVAL(^, D, O). 

Proof. Let G = (U, W, E) be a given connected bipartite graph and label a vertex u £ U. 
Then 

m 

ZXd,o{G) = Y.Za,d,o{KG). 

k=l 

and the values ZA,D,oik, G) can be computed using the EVAL(yl, D, O) oracle by Corollary 

ESI 

The analogous argument labelling a vertex w &W yields the result for Z'^ j^ oiG)- h 

We win show first, that EVAL(A, D, O) < EVAL(^', D', O'). Let G be a given connected 
graph. By equations (|8.6p (j8.7|) and (j8.8|) we have 

Za,d,o{G) = ZXjjq{G) + Z*X^j^q{G) 

= ^A',Di,o'iG)Z^„j~,i,Q„{G) + Z"^, j^i Q,{G)Z*Xii ^Dii ^o"iG) 

By Claim[T]we can compute the values Z^, j^, q,{G) and -^^^j/ o'(.G) using the EVAL(A',D', O') 
oracle. The values Z^„ ^„ q/,{G) and Z^„ ^„ q,{G) can be computed by Claim [1] using the 
fact that EVAL(j4", D", O") is polynomial time computable by the condition of the Lemma. 
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To see that EVAL{A' ,D' ,0') < EVAL{A,D,0) note that by Claim [2 be can compute 

^A,D,oiG) = ^A',D',0'i^)^A",D",0"(.^) 

and 

^A,D,oiG) = ^A',D',0'i^)^A",D",0"i^) 

using an E\AL{A,D,0) oracle. And by Claim [2 using the fact that E\AL{A" ,D",0") 
is polynomial time computable, we can compute Z^„ j~,„ qi,{G) and Z'^,, ^„ q„{G), hence 
^A' D' O'(^) ^^'^ ^A' D' o'(^)' ^^'^ finally 

ZA',D',0'iG) = Z^ij^iQi{G) + ZXi^D',o'{G). 
The proof for non-connected G follows from the above using the fact that 



Za,d,o{G) = \Y^A,D,o{Gi 



i=l 



with Gi, . . . ,Gc being the connected components of G. □ 

8.2 The Proof of Lemma O 

Proof (Proof of Lemma \3.1\) . Let G be a given graph note that if G = iV^E) is not con- 
nected with Gi, . . . , Gk being the components of G then we have 

k c 
i=\ i=l 

This proves (2). To prove (1) note that for hardness we may restrict ourselves to connected 
G. 

Therefore, for some i G [c] fix a component Ai of A and let / C \m\ be the set of 
row/columns indices such that A^ = Ajj. Let G = {V,E) be a connected graph and call 
some vertex z € V the labelled vertex of G. Then by the connectedness of G we have 

ZAAG) = Y,ZA{k,G) 

k£l 
The proof now follows by the Pinning Lemma I8.7[ □ 

8.3 The Proof of Lemma [X^ 

In order to prove Lemma 13.21 it will be convenient to transition from partition functions 
to parity-distinguishing partition functions. How this translation can be performed will 
be described in Lemma 18.151 Once we have determined some conditions on the shape of 
the resulting partition functions the proof of Lemma 13.21 will become straightforward. 



55 



Shape Conditions. Given an evaluation problem EVAL(C, D, O) witli D, O diagonal ma- 
trices of vertex weights and C a connected bipartite matrix with underlying block B. We 
define conditions on the shape of C and D, O. These conditions will be used incrementally, 
that is, we will rely on (C(z + 1)) only if (Cl)-(Ci) are assumed to hold. 

(CI) There are r,m,n £ N, a non-singular r x r-matrix H with entries in { — 1, 1} and 
vectors v £ IRa>0) ^ ^ I^A>o of pairwise distinct entries such that 



B = vnF ® H 



ViWiH . . . ViWnH 
VmWiH . . . VmWnH 



If B satisfies (CI), for convenience, we consider the indices of the entries in B as pairs 

such that -B(^,i),(jyj) = v^w^Hij, for fi £ [m], v £ [n] and i, j £ \r\. We call the submatrices 

v^VpH the tiles of B. 

The diagonal entries of the matrices D and O are vertex weights which by the shape of 

C 

B 

B'^ 



C - \ y 



will be considered with respect to B. As i? is a rm x rn matrix, we group the entries of 
O and D into rm x rm submatrices D^, O^ corresponding to the rows of B and rn x rn 
submatrices D , O corresponding to the columns of B so as to obtain 



D^ \ , ^ / O^ 
Z?^ r^^^ ^ = O^ 



Furthermore, according to the tiles of B the matrix D can be grouped into to m tiles 
£)R,tJ' (^for all fj, £ [m]) each of which is an r x r diagonal matrix. Analogously we group 
the matrix D*^ into n submatrices D^''^ for all v £ [n] and we obtain 

i?«'i ... \ / D^'i ... 

D^= \ ■■ • . . : and Z)'^ = : • . . : 

... L>^'™ / \ . . . L>'^'" 

The matrices O and O are grouped analogously. If B is symmetric then D = D and 
O^ = O^ . We define four more conditions 

(C2) D is a diagonal matrix of positive vertex weights, O^'^ , O^'^ and D + and D — 
are non-negative. 

(C3) The matrix H \s a Hadamard matrix. 



(C4) For all ^ £ [m],u £ [n] there are an a^,a'-^ such that D^'^^ = a^Ir and D 



C,u 



a^Ir. 
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(C5) There are sets A-^, A*^ C [r] such that 

for all /x e [m],zy G [n] there is a (3^,(3^ such that O-^'^ = /3^4;A« and O'^''^ = 

Before we transform a given problem EVAL(A) into the form EVAL(C, D, O) in Lemma 
18. 151 we will exclude some cases from our consideration. That is, we show in the following 
Lemma that EVAL(^) is ^^P-hard unless the block B underlying A satisfies rankabs(-B) = 
1. 

Lemma 8.13. Let A be a symmetric connected bipartite matrix with underlying block B. 
Then at least one of the following outcomes occurs. 

Outcome 1 EVAL(y4) is #P-hard. If B is symmetric, then EVAL(i?) is #P-hard. 

Outcome 2 For some m, n € N there are vectors v E M™ and w G M^ satisfying < wi < 
. . . < Vm and < wi < . . . < Wn and permutations S and n such that 

VmWiS"^^ . . . VmWnS" 

where, for i £ [m] and j G [n], S^^ is a { — 1, l}-matrix of some order rui x rij. If B 
is symmetric then T, = 11. 

Proof. By Lemma l6.4| EVAL(yl) is #P-hard unless rankabs(i3) = 1. Similarly, if B is 
symmetric then EVAL(i?) is ^^P-hard unless rankabs(-B) = 1. 

We conclude that abs(i3) = xy for some non-negative real vectors x,y. If B is sym- 
metric then we can take y = x. To see, suppose x and y are vectors such that xy^ is 
symmetric and let Xi = yt = y/xiin. Note that xiyj = ^/xiinxjyj = xiyj. 

Note that the vectors x and y contain no zero entries. This follows from the fact 
that abs(i?) is a block because B is. Hence if some entry of x satisfies Xj = then 
^i,* = XiX = and therefore B has a decomposition. 

Let V G M™ be the vector of ascendingly ordered distinct entries of x. That is, Vi < Vj 
for all i < j and, for each Xi, there is a j G [m] s.t. Xi = Vj. Similarly, let w be the vector 
of ascendingly ordered distinct entries of y. □ 

Lemma 8.14. Let A be a symmetric n x n matrix of rank r and / C [n] a set of indices 
with \I\ = r. If Aj^: has rank r then the matrix Ajj is non-singular. 

Proof. As ranked/ = rank^ the rows of A with indices in / depend linearly on those 
from /. By symmetry this holds for the columns as well and is still true in Aj. Hence 
rank ^ = rank ^//. □ 

Lemma 8.15. Let A be a symmetric connected bipartite matrix with underlying block Ba 
of rank r. Then at least one of the following outcomes occurs. 
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Outcome 1 EVAL(yl) is #P-hard. If Ba is symmetric, then EV'AL{Ba) is #P-hard. 

Outcome 2 There is a connected bipartite matrix C, whose underlying block B is size 
mr X nr for some m and n, and diagonal matrices D and O which satisfy conditions 
(CI) and (C2), such that 

EVAL(C, D, O) = EVAL(y4). 

The matrices D and O consist of mr x mr submatrices D^, O^ and nr x nr subma- 
trices D ,0 such that 

D=i j^c ) andO=[^ ^ ^c 

C,D and O can be computed in time polynomial in the size of A. If Ba is symmetric 
then so is B. Also D^ = D^ , O^ = O^ and 

EVAL(B, D^, O^) = EVAL(Ba). 

Proof. Suppose that the matrix A does not give Outcome 1 in Lemma 18.131 Let S and 
n be the permutations from Lemma 18.131 and let $ be the permutation on the rows 
of A that apphes S to the rows of Ba and apphes 11 to the columns. Let A = ^4$ $. 
Note that EVAL(A) = EVAL(2). Also, the block underlying A is {Ba)y.,ii-> which we 
denote B. Note that B is symmetric if Ba is symmetric, since S = 11 in that case and 
EVAL(Ba) = EVAL(5). By Lemma EH] there are m,n G N such that 

ViWiS^^ . . . ViWnS^^ 

B 

for vectors v G M™, w G M^ of positive pairwise distinct reals and { — 1, l}-matrices S of 
order m^ x nx- Let 

5^1 ... 5I" 

S 

oml or, 

For convenience, we consider the indices of the entries in B as pairs such that Bu^i)ix,j) — 
Vf^w\S1^, for (k, A) G [m\ x [n] and (i, j) G [m^] x ['^a]- Entries and submatrices of S will 
be treated in the same way. 

First we shall see that we may assume that every pair of rows (or columns) of S is either 
orthogonal, or they are (possibly negated) copies of each other. 

Claim 1. Outcome 1 occurs unless for all k, A G [m] and i G [?n.K],i G [m\] 

either {Sf^^, S^^ = for every u G [n] 
or there is a s G { — 1, +1} such that S'^'^ = sSj''^ for every z^ G [n]. 
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The analogues holds for the columns of S: for all k, A G [n] and i £ [t^k], J G ['^a] 

either (5^^, S'^ ■) =0 for every /i G [m] 
or there is a s G {—1, +1} such that S'^'^ = sS^ ■ for every fi G [m]. 



(8.10) 



Proof. Let p G N be odd. By p-thickening and subsequent 2-stretching we obtain a reduc- 
tion 

EVAL(^') < EVAL(I) 

for a matrix A' = (1(p))2 which contains submatrices B'^p\B^p>Y and {B^p'^)'^B^p\ The 

same reduction gives EVAL((i?(^') ) < EVAL(S) if B is symmetric. We will give the proof 
of equation (j8.9p by focusing on B^P' [B^'^')'^ . The analogous argument on [B^^'Y' B^p> 
yields equation (jS.lOp . 

Let C = B^P\B'^P^y . For k, A G [m] and i G [m«;], j G [mx] we have: 

(v,k) v=l 

Note that by 2-thickening we have a reduction EVAL(74") < EVAL(A) for a matrix A" = 
{A')^'^\ This also gives a reduction EVAL(C'(2)) < EVAL(5) if B i^ symmetric. The 
matrix A" has only non-negative entries and contains the submatrix C'^^ The result of 
Bulatov and Grohe [6j implies that EVAL(C'(2)) and EVAL(A") are #P-hard, in which 
case, Outcome 1 occurs, if C^"^' contains a block of row rank at least 2. We shall determine 
the conditions under which this is not the case. 

A 2 X 2 principal submatrix of C^^' , defined by {k^i), (A,j) has determinant 



det(«;,i),(Aj) 



'^(A,i),('^,i) ^(A,i),(A,i) 

We have 

/n X2 

\K,i),(K,,i) 



{C{^,i),(K,i)C(x,j),{x,j)f - iC(^,i),(x,j)Y 



^£),(M = <'{E^l'is^,:,s^p] =vtp{Y,wlPn. 



\u=l / \i^=l 

and therefore 

det(.,),(A,,) = vt^vf lif^y^l'r^u] - f E ^^^(5.':*^ SJ;:) 

This determinant is zero iff there is an s G {—1,1} such that {Sf^,S^'^^) = suu for all 
V G [n\ which implies Sf^ = sSj'^ for all z^ G [n]. By equation (jS.lip and Lemma 16.51 we 
further have C}-^ (A "i ~ ^ ^°^ arbitrarily large p iff {S'^';^, S^'^) = for all z^ G [n]. ^ 
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Assume from now on that Equations (j8.9p and (jS.lOp hold. The next claim states that 
the rank of each tile of S equals the rank of S itself (which is equal to r, the rank of B 
wich is the rank of B^). 

Claim 2. rankS" = rank 5'^'^ for all (k, A) G [m] x [n\. 

Proof. Equation (j8.9p implies that rankS""^ = rankS"''' for all k S [m\ and /u,z/ G [n]. 
Combining this with equation (jS.lOp we obtain rankS'''^ = rankS'^'' for all k, A G [m] and 
IX, u & [n]. 

Therefore it suffices to show that r = rankS = rankS*^^ holds. Let S"*^ denote the 
matrix 

/ 511 

\ ami 

Let / be a set of row indices with |/| = rank S = r such that the set {Sj^* | i G /} is linearly 
independent. By equation (18. 9p we have (S'j,*, S'j,*) = and (5*^, S**!^) = for all i / j G /. 
Hence, S*^ has rank r. As S^^ is a mi x ni matrix there is a set J C [ni] s.t. the columns 
of S*^ with indices in J form a rank r set. Equations (jS.lOp implies {SW, 5'^l) = for all 
i ^ j £ J. This proves the claim. -I 

Claim [2] has strong implications on 5" ( and B). It implies that for all (k,A) G [m] x [n] 
there are sets K(^^^X)^^{k,X) of cardinality r such that S'^ ^ is non-singular. By 
equation (j8.9p we take, without loss of generality, K^^i^x) = ^{k,a') foi^ ^^ ^ ^ ["^] ^^id 
A, A' G [n]. Analogously, equation (jS.lOp implies -^(k,a) = -^(k',a) for all k, k' G [m] and 
A G [n]. Therefore, there are sets of indices Ki, . . . , Km and Li, . . . , L„ each of cardinality 
r such that the matrix 

SrIlx is non-singular for all (k. A) G [m] x [n]. (8.12) 

If i?A is symmetric then B is symmetric and we may assume, by Lemma 18.141 that 
i^K = -^K for all ^ S [m]. But there is more we can infer from Claim [H namely the above 
non-singular subtiles of each tile are (up to row-column negations and permutations) equal: 



Claim 3. For all k G [m] and A G [n] the sets K^^ and L\ have orderings 

Kf, = {K,i, . . .,k^,r} and Lx = {ix,i, ■ • .,ix,r} 

and there are families {r^ : [r] -^ {-1, l}}Ke[„^] and {rf : [r] -^ {-1, l}}Ae[n] of mappings 
such that: 

Slli,^, = r^{a)r^{b)S^kt,i..^ for all (k. A) G [m] x [n], a, 6 G [r]. 

If B is symmetric then S is symmetric and i^^ = ^k and r^ = r^ for all k G [m] . 
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Proof. As Sj^ ^ = SkiLi and rank 5^^ = rank 5, equation (j8.9p implies that every row in 
S is either a copy or a negated copy of a row in Sk^*- Fix an arbitrary ordering Ki = 
{ki^i, . . . , ki^r}- As 5^\ ^ has rank r for all k € [m] there is an ordering {/ck,!, • • • , /c^.r} 
and, for every a G [r], an Sa G {-1,+!} such that 5*^^^,, = ■Sa<S'^^^_*- Let r^(a) = Sa- 
Then S^-^ ^ = T^{a)S^^ ^ for all a G [r]. Equation ()8.9p implies that this extends to 

Sit,. = ^^i^)S^t,* for ah a G [r],K G [m], A G [n]. 

An analogous argument on the columns of S using equation (jS.lOp . yields orderings of 
the sets Lx and mappings tx such that 



S:},,, = r^ib)S:l^, for ah b G [r],K€ [m], A G [n]. 
Combining both finishes the proof of Claim [3j H 

For K G [m], let vr^ be a permutation of [m^] which satisfies 7r^(o) = A;K,a for a G [r]. 
For A G [n], let vr^ be a permutation of [nx] which satisfies vr^(a) = £x,a for all a (^ [r]. 

Let 5"^ be the result of the permutations irf} and vr^ when applied to S that is 
5"''* := {S'^^)^R ^c . Let -B be the matrix defined by Bi^^^iy^x.j) = VnWxS^j and let S be the 
matrix defined by S'(K,i).(A.j) = '5'f^- Let A be the bipartite matrix with underlying block 
B. Note that EVAL(i) = EVAL(2) = EVAL(A). The definition of these permutations 
implies that B is symmetric if B is symmetric (which is true if B^ is symmetric). In this 
case, EVAL(S) = EVAL(i3) = EYAL{Ba). Equation (|8l^ simphfies to 

S'f^.ir^i is non-singular for all(K;, A) G [m] x [n] (8.13) 

and Claim [3] implies furthermore that 

5^,1 = r,^(a)rf (6)5f, for all (k. A) G [m] x[n],a,be [r]. (8.14) 

We consider the twin-relation on A now. As A is bipartite, the equivalence classes of 
this relation induce collections of equivalence classes separately for the rows and columns 
of B. Furthermore, as -B(^ j) (;^j) = v,^wxS1^ and the values Vi are pairwise distinct and 
positive, two rows corresponding to different vi values are not twins. This is similarly true 
for the columns of B. Hence, the equivalence classes of rows can be grouped into collections 
Xi, . . . ,Tm and the equivalence clases of columns can be grouped into collections Ji, ■ ■ ■ ,Jn 
such that, for every k G [m], the collection X^ contains the equivalence classes of rows in 
the submatrix 

T^* := ( v^wiS^^ . . . v^WnS^"" ) 

of B. By equation ()8.9p and equation (I8.13P every row in T^* is either a copy or a negated 
copy of a row in (T''*)r^i^. Moreover, every two i ^ j £ [r] belong to different equivalence 
classes by equation (I8.13p . 

We may therefore assume, without loss of generality, that the collection X^ consists of 
classes Pf *, . . . , Pp'* and iVf*, . . . , A''^^* such that i G Pj^* for ah i G [r]. Furthermore, the 
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sets N!j^* account for the possible negated copies of rows in (T'^*)rrU and therefore some 
of these sets may be empty. But for all i G [r] if A'^^'^* is non-empty then all a G N^* are 
indices of negated copies of rows from P^* . 

Similarly, the collection Jx of equivalence classes of columns corresponds to the subma- 
trix 

'T^*A I 

V VmWxS"'^ 

of B. By equation (|8.1Up every column in T*^ is either a copy or a negated copy of a 
column in {T* )m- Moreover, by equation (j8.13p every two i ^ j £ [r] belong to different 
equivalence classes of the twin relation. 

We may assume that the collection J'x consists of classes Pi'^, ■ ■ ■ , P*^ and A^-j*''*, . . . , N*^ 
such that i G P*^ for all i G [r\. The sets N*^ account for the possible negated copies of 
columns in (r*^)^r^i and therefore some of these sets may be empty. But for all i G [r] if 
N* is non-empty then all a G A^* are indices of negated copies of columns from P* ■ 

Note that if B is symmetric the above definitions directly imply that m = n and, for all 
/i G H, T^ = J"^. Also, we can take P^"* = P-^ and N^* = N*'' for all i e [r]. 

Application of the Extended Twin Reduction Lemma fS. 31 according to these equivalency 
classes therefore yields an evaluation problem EVAL{C,D,d) = EVAL{A){= EVAL(y4)) 
such that the block B' underlying C satisfies 

( viWiSl^][r] ••• ^i^n^HH ^ 
B' = 



\ VmWlbr-ir-: . . . VmWfiJl^U^] I 



That is, B' is an mr x nr matrix and D and O are diagonal matrices of vertex weights of 
order mr + nr. Grouping these vertex weights according to the rows and columns of B' to 
which they correspond, we obtain 

for mr X m,r diagonal matrices D^,0^ and nr x nr diagonal matrices D^ ,0'~' . Their 
structure corresponding to the tiles of B' in turn is 

i?«'i . . . \ / Z^^'i ... 

D^ =\ ■■ • . . : and Z)'^ = : • . . : 

... L>^'™ / \ . . . D^'" 

which holds analogously for O such that the D^'^, O^'^, D*^'*^, O^''^ for all /i G [m], v G [n] 
are r ^ r diagonal matrices. The definition of these matrices according to the Extended 
Twin Reduction Lemma 18.31 is then, for all ^ G [m], z^ G [n], i,j G [r], given by 

(8.15) 





K*\-\Nr\ 


and Dfj" = 
and df^f = 


\P*'^\ _|_ 1 JV*" 

l-Tj 1 ^ I^^J 

l-^J 1 I^^J 
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If Ba is symmetric then B is symmetric and D = D . Also, B' is also symmetric and 
EVAL(^',L>^,0^) = EVAL(S). 

Clearly, the matrix D is a diagonal matrix of vertex weights whose diagonal is positive 
as the sets P^* and P*^ are non-empty by definition for all k G [m], A G [n] and i G [r]. 

By Equation (|8.14p . for all (k. A) € [m] x [n], the matrix 5'tu^i is - up to negations of 

rows and columns - just a copy of the matrix S}\,y However, the diagonal entries of O 
given by equation (j8.15p may be negative in some cases. To satisfy condition (C2) we 
therefore define mappings /> : [r] — > { — 1, 1} and 7 : [r] — > { — 1, 1} by 

P{^ = [-' '^'^5^<° and ,{J) = [-' ^^^^<' 

I 1 , otherwise I 1 , otherwise 

We will use these mappings below to "transfer" the signs of diagonal entries of O^'^ and 
O '^ to B'. Note that /> = 7 if Ba is symmetric since O = O in this case. Define 
matrices 5'f^u^i by applying row and column negations according to these mappings, that 
is 

S^H = p{aUb)T^{a)T^{h)S^al for ah (k. A) G [m] x [n], a, & G [r]. (8.16) 

By equation (I8.14p . we have the following for all (k, A) G [m] x [n] and a,b ^ [r]: 

Thus 

Tf{a)T^{b)S'^, = p{aMb)Tf{a)T^{b)Tf{a)r^{b)S'^, = 5,-,i 

But, by their definition in Claim[3l the mappings r/^ and rf satisfy t/^(«) = t'~" {i) = 1 for 
all i G [r]. So the above equation gives 5^^^ = S'l\ for all (k, A) G [m] x [n] and a,6 G [r] 
so Srir^,] = Sy^,, for all (k, A) G [m] x [n]. Define H := S}},,. Let B be the matrix defined 
by ^(K,i),(Aj) = Vf^wxHij so 

B{K,i),{x,j) = v^wxSij 

= v^wxp{t)l{3)r^ii)T^{j)Sf;; 
= p(i)7(j)r^(i)rf (j)B(^,i)_(Aj). 

Let C be the symmetric bipartite matrix with underlying block B. For k G [m], A G [n] 

and i, j G [r], let 



^ff = Pii)r^(^O^t and og^ = 7(j>A^(i)0£'. 



Let O^ be the diagonal matrix with tiles O^''^ for k G [m] and O*^ be the diagonal matrix 
with tiles O ' for AG \n]. Let O be the matrix 



O^ 
O^ 



O = . . ^c 
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Since (as noted above) T{^(i) = t (i) = 1 for all i G [r], the matrices O '^ and O '^ are 
non- negative. 

The Row-Column Negation Lemma 18.41 implies 

EVAL(C', D, 6) = EVAL(C7, D, O). 

The block B satisfies (CI) and the matrices D and O satisfy (C2). The definitiond of 
D and O in equation (jS.lSp and the definition of O implies that D + O and D — O are 
non-negative as required. If Ba is symmetric then B and S are symmetric so t^ = t^ 
and vr^ = vr^ so S and B are symmetric. Since p = j , B is also symmetric. So the 
Row-Column Negation Lemma 18.41 implies 

EYAL{B', D^, 6^) = EVAL(S, D^ , O^). 

Furthermore it is easy to see that all operations performed to form C, D, O from the 
matrix A are polynomial time computable. This finishes the proof. □ 

The remainder of this section relies on a gadget which consists of arrangements of paths 
of length 2. These paths affect the matrices C,D,0 in a similar way to 2-stretching. 
It is therefore convenient to have a look at the effect this operation has. Clearly 2- 
stretching yields EVAL{CDC,D,0) < EYAL{C,D,0). If B is symmetric it also yields 
E\AL{BD^B,D^,0^) < EYAL{B , D^ , O^) . 

Assume that C and D, O satisfy conditions (CI) and (C2). Recall that B = vnF (8> H 
holds for the block B underlying C . Furthermore the matrix CDC contains the subma- 
trices BD^B'^ and B'^D^B and 



Wlfl 



HiE:=i^iD''nH^ ••• v,v^Hiz:=,wiD^nH 



T 



bd^b'^ 



with analogous analysis of B^D^B we have 

(n \ j m 

H(^ wId'^''')H^ j and B^D^B = wuF ® H'^C^ vId^''')H 
v=\ ) \ ^=1 

(8.17) 

We define a reduction template T{t,p,q) which will be used in the proofs of Lemmas 
18.171 and 18.191 Let P{t,p) be a graph constructed as follows. Start with an edge with a 
distinguished endpoint a. Then perform in succession a f-thickening, then a two stretch, 
and finally a p-thickening. (Informally, there is a vertex b connected to o by t many length 
2 paths such that all edges in those paths have multiplicity p.) 

The reduction T{t,p,q) works as follows. In a given graph G = {V,E), we 2-stretch 
each edge e G E and call the middle vertex Ve- We attach q disjoint copies of P{t,p) by 
identifying their terminal vertices with Ve. Figure [3] illustrates the construction. 

Recall that M o N denotes the Hadamard product of matrices M and N . 
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U V U V 

Figure 3: The gadget templates r(l, 3, 2) and T(2, 2, 1) 



Lemma 8.16. Suppose C andD,0 satisfy (CI) and (C2). At least one of the following 
outcomes occurs. 

Outcome 1 EYAL{C,D,0) is i^P-hard. If B is symmetric, then EVAL{B,D^,0^) is 

#P-hard. 

Outcome 2 Fort,p,q G N andp' = 2p+l and q' = 2q there are rxr matrices @ = @{t,p') 
and S = E{t,p') defined by 

{Ip') l^% S QR,f, if t is odd 



(7,^)'-E- 



u=l 



ipi ( D'-'''^ , if t is even 
" ■ 1 0^> , ift is odd 



for positive constants 7I; and 7^, depending on p' . 

The reduction T{t,p',q') yields E\AL{C AC , D , O) < EYAL{C, D,0) for a diagonal 
matrix 



A^ 
A^ 



A = A{t,p',q') = ( „ ^r' 1 cind a matrix CAC 



BA^B'^ 

B^A^B 



A is a diagonal rmxrm matrix of rxr tiles A '^ = Vff '^ D '^oQ'^'i > for all /x G [m]. 
A^ is a diagonal rn x rn matrix of rxr tiles AP'^ = Wy '^ D*^'^ o'E^'^ > for all u ^ [n]. 
If B is symmetric then the same reduction yields 



. R T3 nR nR 



R r)R\ 



EVAL(SA"5,L>^,0") < EYAL{B, D'^,0 



Proof. Let p',q' be as above. 

Claim 1. Either Outcome 1 occurs or there are constants 'yS and 7^ depending on p' such 
that 

B(p')d''{b'^p'Y = (vv^)^"'^ l3lr and (i?(p'))^Z)^i3(p') = {ww^)(p'^ 7S/, (8.18) 
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Proof. We have EYAL{C^p'^DC^p'\D,0) < EVAL(C,D,0) by p'-thickening followed by 
2-stretching. If 5 is symmetric this also yields EVAL(5(p')D^S(p'), L>^, O^) < EYAL{B, D^,0^). 
The matrix C7(P')Z)C7(P') contains submatrices X := B^P'^ D^ {B^P'^f andY := {B^P'^fD^B^P'\ 
We show the first part of equation (jS.lSp by an argument based on the matrix X. The 
second part then follows analogously using Y. (Recall from (CI) that D = D when B 
is symmetric, in which case X = Y.) 

Define U = J2u=i wlD^'". By equation ([gTTl) we have X = {vv'^)p' ® {HUH'^) . There- 
fore, if ahs{HIlH^) contains a block of row rank at least two then X does. 

As iif is a {—1, l}-matrix we have {HIlH'^)ii = tr(n) for all i E [r] and the trace of 11 
is positive. Furthermore \{HIlH'^)ij\ < tr(n) for all j ^ i hy the non-singularity of H. 
Hence, we obtain a block of rank at least 2 in ahs{HIlH^), if there is a non-zero entry 
{HIlH^)ij for some i ^ j G [r]. The proof follows with 7^ = tr(n). -| 

For convenience, let T = T(t) denote the matrix D, if t is even, and O otherwise. 

Recall the reduction template, let (fi, i), {k, k) £ [m + n]x [r] denote the spins of Ve and 
b. 

The diagonal (/i, i) entries of A correspond to the partition function of the reduction 
template with vertex Ve fixed to {n,i)- Therefore, for fi G [m] 

(m r \ *? 

K=l fc = l / 

(■m r \ 1 

K=l fe = l / 

(m \ a' 

where the last equation follows from Claim 1. Similarly, for i^ G [n] 

/ n r \ <?' 

A^ = D?fiEET^:^'(cip'^Dc('X.u.^)] 

\K=lk = l / 

/ n r \ <?' 

= DfflY^Y.T^f{iB^p'rD^Bipx, 

\k=1 k=l 

'Dffih^^YtwiP'Tl 



i),{K,k) 
K=l k=l / 



n \1 



With and H defined as in the statement of the Lemma the proof follows. □ 

Lemma 8.17. Let C and D,0 satisfy (CI) and (C2). At least one of the following 
outcomes occurs. 
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Outcome 1 EYAL{C, D,0) is #P-hard. If B is symmetric, then EYAL{B, D^,0^) is 

#P-hard. 

Outcome 2 Conditions (C3) and (C4) are satisfied. 

Proof. The ^^P-hardness part will be shown using a gadget construction T{2,p',q') with 
p' = 2p+l and q' = 2q for p, g G N. By Lemma[836]this yields a reduction EVAL(CAC, D, O) < 
EVAL(C,Z?,0) such that CAC contains submatrices BA'^B'^ and B'^A^B. If B is 
symmetric then EVAL{B A^ B , D^ , O^) < EYAL{B , D^ , O^) Focusing on B^A^B we 
will prove (C3) and the part of (C4) which claims that D '^ = a]]^Ir- The proof for 
D '"^ = a^ Ir then follows by analogous arguments based on BAB. 

Recall that by the proof of equation (j8.17p we have B^A^B = {wnF) (H'^A'H) for 
an r X r diagonal matrix A' defined by 

m I m \ 

^=1 \^i=l I 



with eW = @{2,p',q') = (7^)2 -^vlP' • D^'^'. (8.19) 

If ahs{H^ A' H) contains a block of rank at least 2 then abs(i?A^i?) does. So, if 
ahs{H'^ A' H) contains a block of rank at least 2, then Outcome 1 occurs by Lemma |6.4[ 

By the definition of H A'H, we have {H A' H)i^i = tr(A') for all i G [r] and this trace is 

positive by the definition of A'. Therefore, every principal 2x2 submatrix of abs{H^ A'H) 

has the form 

tr(A') \{H'^A'H)i,j 

{H^A'H)j,,\ tr(A') 

As H is non-singular \[H A'H)ij\ < tr(A') for all i ^ j £ [r] and therefore, every such 
submatrix has non-zero determinant. Furthermore, such a submatrix is part of a block if 
(H A'H)ij 7^ 0. Therefore we have Outcome 1 if we can show that {H A'H)ij ^ for 
some i j^ j £ [r] and some p,q £ 'N. 

Assume therefore that {H A'H)ij = for all i ^ j £ [r] and all p, g G N. The remainder 
of the proof is to show that in this case conditions (C3) and (C4) are satisfied. 

Let i?p,g,i = ^2"^^! vfF''^'^'^Df;f for all i £ [r]. Note that A^ = t?p,g,i0W -'i. We define an 
equivalence relation ~ on [r] by letting i ~ j if and only if D^ f^ = D- '-^ for all fi £ [m]. 
Let I be the set of equivalence classes. We will use the notation Dj '^ to denote the value 
DJ^f for i E /. 

Recall that the values v^ in the definition of ??p,q,i are pairwise distinct and non-negative. 
Lemma 16.51 implies the following, for all i,j £ [r]: 



(For all p and q^-dp^q^i = ^p,q,j) iff i ~ j. (8.20) 
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(For all p,eWi,i = e^jj) iff i ~ j. (8.21) 

We use the notation 'i?p,q,/ to denote the value 'i?p,g,j for i £ I. Similarly, we use the 
notation Qpj to denote the value Of- for i e I. 

For i,j £ [r] define sets Vij = {k £ [r] \ H^iH^j > 0} and Mij = {k £ [r] \ HkiH^j < 
0}. 

Then we have 



{H A'H)ij = 2_^Hk^iHkjA'i^f^ = 2_^Hk,iHk,j^p,q,k[&^^^ 



,9 

'k,k 

fc=l fc=l 



E ^p,9,fc(0'"'fe,fc)'^' - E ^p,^A&^'''^i,iy' ■ 



Then 



E^p!/ E ^p.^-^- E ^p.'/-^ 

/ei \keinVij leinMi^ 

E^p.'/-^0p!/(l^nP„-|-|/nAA,,,|) 



Claim 1. Suppose that (H^A'H)ij = for all i ^ j £ [r] and all p,q £'N. Then there is 
a J G I such that \J n Vijl = \J nMij\ for all i ^ j £ [r]. 

Proof. Choose p sufficiently large that there is a unique J £ 1 maximising Yl^=i '^J' ^/^ 
For this p, and for any q £'H and / G I, we have < ^p,qj < '&p,q,j and < Qp^i < @pj- 

Now consider i ^ j £ [r]. For ah / G I, let cj = \I D Vijl - \I D Mij\. 

Since {H'^A'H)ij = 0, for all g G N, 

= E^^^p.^-^0?/ 

lei 

= Cj1^p,q,J&llj + E ^I'^P,<l,I®l'!l 
/GI\{J} 

= C.+ E c,^(|^)". (8.22) 



As g tends to infinity, the sum tends to so cj = 0. 
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Assume now that {H A'H)ij = for all i ^ j ^ [r] and p,q G N. Fix J G I such 
that I J n Vij\ = I J n Mij\ for all i ^ j £ [r]. Recall that Hj:,, denotes the the submatrix 
of H consisting of the rows of H with indices in J. For each pair i ^ j G [r], the fact 
that \JriVij\ = \JnMij\ implies ((-ffj,*)*,i, (-f^j,*)*,j) = 0. Hence, the columns in Hj,, are 
pairwise orthogonal. Since the rank of H is r, this implies that \J\ = r. Now since the 
rows of i/^ are pairwise orthogonal, we have H^H = rl^ so the inverse of H^ is r~^H. 
As right inverses of matrices are also left inverses, we have r~^HH = Ij. and therefore H 
is a Hadamard matrix and we have proved condition (C3). 

Finally, J = [r] implies that D^l^ = D- '-^ for all i,j G [r\. Equivalently, D '^^ = a^Ir 
for some appropriate a^. This proves (C4). □ 

We call a diagonal matrix D pre-uniform if there is a non-negative d such that all 
diagonal entries Dii of D satisfy Di^i G {0, d}. An important technical tool in the last step 
of our proof of conditions (C1)-(C5) will be the following Lemma. 

Lemma 8.18 (Pre-Uniform Diagonal Lemma). Let H be a non- singular rxr {—1, 1}- 
matrix and D he an r x r diagonal matrix with non-negative entries in M. // D is not 
pre-uniform, then there is a p ^'H such that ahs{H D^p> H^) contains a block of row rank 
at least 2. 

Proof. Note that, if the diagonal of D is constantly zero then D is pre-uniform. Assume 
therefore that there is some positive diagonal entry in D. Define B := HD^^'H^, K := 
{k G [r] I Dk^k > 0} and s := \K\. Hence, for i,j G [r], 

(8.23) 
D 
That is, for every / C [r], we have Bjj = Hj^kD^ j^{Hj^k)'^ ■ Fix a set / C [r] such that 
|/| = s and the matrix Hjk has rank s. Since Hjk is non-singular, every 2x2 principal 
submatrix of Bjj has non-zero determinant. To see this, note that, by equation (18.230 we 
have Bi^i = ti{D^^) for all z G / and this trace is positive. Then every such principal 
2x2 submatrix has determinant 

tr(Z)g^) \{Hj,KD^p^{Hj,Kf),, 

\iHj,KD^P,,{Hj,KVU tr(Z?g^) 

and by the non-singularity of Hjk we have \{Hj^kD^ j^{Hj^K)'^)i,j\ < ^^{^k k) (compare 
equation (18.230 ). Hence the above determinant is non-zero. 

Assume that, for all p G N, there are no non-trivial blocks in Bjj, i.e. Bij = for all 
i ^ j £ I. We will show that this implies that D is pre-uniform. 

Fori,j G / define the sets "Pj J := {k G K \ Hi^kHj,k = 1} andAAjj := {k G K \ Hi^kHj^k = — !}• 
That is, Vij and Mij form a partition of K. Therefore, for i,j £ I we have 



^ij = Yl ^i,kHj,kDl^. = Y^ D'li^ - Y^ Dl 



k,k- 



fc=i keVi^j k&Ni,j 
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Partition K into equivalence classes J such that i,j £ K are in the same equivalence 
class iff Di^i = Djj. Let J be the set of these equivalence classes and for each J £ J 
define Dj := Dj^j for some j S J. We have 






Jf- 



As the Dj are positive and pairwise distinct Lemma 16.51 implies that, (for all p we have 
Bi^j = 0) iff (|JnPjj| = iJnA/ijI for all J). By our assumption that this is true for 
all i ^ j €z I we see that the s x | J| matrix Hjj is orthogonal which implies \J\ = s. In 
particular, J = K and Dk,k 

Lemma 8.19. Let C and D,0 satisfy conditions (CI) - (C4). At least one of the fol- 
lowing outcomes occurs. 

Outcome 1 EVAL{C,D,0) is #P-/iard If B is symmetric, then EVAL{B,D^,0^) is 

#P-hard. 

Outcome 2 Condition (C5) is satisfied. 

Proof. We will use reduction template T{l,p',q') with p' = 2p + 1 and q' = 2q for p,q £ 
N. By Lemma KTE\ this yields a reduction EVAL(CAC, D, O) < EVAL(C, D, O) such 
that C/S.C contains submatrices BA'-^B^ and B^A^B. If B is symmetric then it yields 
the reduction EYAL{B A^ B , D^ , O^) < EYAL{B , D^ , O^) . We base our argument on 
B^A^B to prove that O^'f" = l3^Ir-xR for all n G [m] and some (5^ and A^ C [r]. The 
analogous argument on BAB then yields the result for the submatrices of O . 

Recah that by equation (j8.17p we have B^A^B = {wuF) ® (H^A'H) for an r x r 
diagonal matrix A'. With 

m 

eW = e{i,p') = ip'-Y^ yp^ ■ o^-'^ (8.24) 

the r X r diagonal matrix A' is defined by 

m j m \ I m \ 

The last equality holds by condition (C4). Taking 

m 

^'■=Y1 <''''^^«M we have A' = ^G^^'^'^ (8.25) 

If ahs{H A'H) contains a block of rank at least 2 then abs(i? A S) does. So, if 
ahs{H^ A' H) contains a block of rank at least 2, then Outcome 1 occurs by Lemma |6.4[ 
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By the definition of H^A'H, we have {H A'H)i^i = tr(A') for all i G [r] and this trace 
is non-negative by the definition of A'. Therefore, every principal 2x2 subniatrix of 
ahs{H^A'H) has the form 



tr(A') 
\{H^A'H] 



j,*i 



\{H'A'H)i,,\ 
tr(A') 



As H is non-singular \{H A'H)ij\ < tr(A') for all i ^ j € [r] and therefore, every such 
subniatrix has non-zero determinant, if tr(A') is positive. Furthermore, such a submatrix 
is part of a block if {H'^ A'H)ij ^ and tr(A') ^ 0. Therefore we have Outcome 1 if we 
can show that {H A'H)ij ^ and tr(A') ^ for some i ^ j G [r] and some p, g G N. 

Assume therefore that either {H'^A'H)ij = or tr(A') = for all i ^ j G [r] and 
all p,q S N. The remainder of the proof is to show that in this case condition (C5) is 
satisfied. 

Recall that by equation (j8.25p the value 'd is positive for all p,q £ 'N. Therefore A^ ^ = 
iff eWi,i = 0. 

Claim 1. There is a po ^ N such that for all p > po and all i £ [r] we have 

G^i.i = iff {off = for all /i G [m]). 



Proof. For each i G [r] application of Lemma 16.61 to equation (j8.24p yields that there is a 
Pi such that for all p > Pi we have 

G^i.i = iff {off = for all p£[m]). 



The Claim follows with pg ■= maxjpi, . . . ,pr}. 



1(2) 



Claim 2. Let p £ N. If {H'^A'H)ij = for all i / j G [r] and all g G N then 6^^ ^ is 
pre- uniform. 

Proof. Define n = (0W)(^). Then all entries of LI are non-negative and H^'^) = (0[p1)('? ). 

With H^A'H = ^{H'^Q^P^^'^'^H) = ^{H'^Ii^'i'i H) the claim follows by the Pre-Uniform 
Diagonal Lemma 18.181 H 

Claim 3. There is a p= G N such that for all p > p= and all i,j£ [r] we have 

&% = ©'"''L iff (Off = off for all p£[m]). 

Proof. The backward direction holds for all p G N. Fix i,j G [r]. By equation (|8.24p the 
equality Q'-^'i^i = © jj implies 












,3 
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By Lemma 16.61 we either have a pij such that for all p > pij 

el^l ■ . = eW',. iff Oi;t = off for all /x G [. 
or there is a p'^- such that for all p > pj- 

0'"'m = ©'"'L- iff Og'^ = -off for all /x G [ 



m\. 

However the second possibility would particularly imply that O^l < for some i € [r] 
which was precluded by condition (C2). Therefore the first possibility holds with pij. 
The Claim now follows with p= = maxjpjj- | i, j E [r]}. ^ 

These claims now enable us to finish the proof. Note first that condition (C5) is satisfied 
with A^ = if off = for all n G [m] and i£ [r]. 

Assume therefore that O^ has non-zero diagonal entries. Fix values po,P= S N according 
to Claims [1] and [3] and define p = max{po5P=}- This implies tr(A') ^ 0. To see this note 
that there is some i E [r] and some fi E [m] such that Oj ■ ^ which by our choice of p 
implies 9^^,^ / 0. 

By our assumption, tr(A') ^ implies {H^A'H)ij = for all i ^ j ^ [r] and all g E N 

which by Claim [2] yields pre-uniformity of G^^^ 

Define A^ := {i E [r] \ 9^^,^ / O}. By the pre-uniformity of 9^^^^ Claim [3] imphes 
that, for each fi E [m] and every i E A there is a Pf such that O^ f = Pf. Furthermore, 

Claim [1] implies that for each /i E [m] and every z E [r]\ A we have O^ f = 0. This 
finishes the proof. □ 

8.3.1 Putting everything together. 

We are now able to prove Lemma 13.21 

Proof (Proof of Lemma \3.2\) . Bipartite A. Consider first the case in which A is bipartite. 
By Lemmas 18.151 18.171 and Lemma 18.191 the evaluation problem EVAL(74) is #P-hard 
unless EVAL(^) = EVAL(C,D,0) for matrices C,D,0 satisfying conditions (C1)-(C5). 
C is a symmetric bipartite matrix with underlying block B. Conditions (C1)-(C5) 
imply that B = vuF ^ H, D^ = D^" Ir, D^ = D^" /^ and O^ = O^" It-kr-, 
O = O (8) Ir;k'^ fo'^ diagonal m x m matrices D and O defined by Df,j^ = af 
and of = Pf for all jj, E [m]. The n x n diagonal matrices D'-' and O'-' are defined 
analogously in terms of a^ and P^ . Then for 



O" = ( """ ° ■^ a,.d O" = ( O"" » ■^ and C" ^ ° 



vw 



T 



j^c" J --- " - V O^" J "^ ~\wv^ 

the problem EVAL(C",I?", O") is polynomial time computable by Corollarv 18. 2[ Note 
that D + O and D — O are non-negative by condition (C2). Hence with M, A being the 
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bipartisation of F,A^ and A*^ we have EVAL(C, D, O) = EVAL(M, /2r,/2r-;A) by Lemma 

Em 

Non-Bipartite A Now suppose that A is not bipartite. Let M be the bipartisation of 
A. Recah that this is a matrix of the form 

By Lemmas I8.15( 18.171 and Lemma 18.191 the evaluation problem EVAL(A) is #P-hard 
unless there are matrices C,D,0 with block B underlying C satisfying conditions (Cl)- 
(C5) such that EVAL(^) = EVAL(5, D^, O^) 

Conditions (C1)-(C5) imply that B = vv'^^H, D^ = D^"^Ir, and O^ = O^" ^I^.j^r 
for diagonal m x m matrices D^ and O^ defined by -D^ = a^ and Oj}^ = (5^ for all 
^ G [m]. 

Then the problem 'ENKL{vv ,D ,0 ) is polynomial time computable by Corollary 
MM Hence we have EYAL{B , D^ , O^) = EVAL(^,/^,I^.^fl) by Corollary ETH 

Finishing the Proof: It remains to state the polynomial time computability. Note 
that conditions (C2) (C5) are straightforwardly checkable in polynomial time and for 
(CI) this follows from Lemma l8.15i □ 
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